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(g) Method and apparatus for multiple access between transceivers . in wireless communication using 
OFDM spread spectrum. 

(g) A method for allowing a number of wireless transceivers to exchange information (data, voice or 
video) with each other. A first frame of Infomnation is multiplexed over a number of wideband frequency 
bands at a first transceiver, and the Information transmitted to a second transceiver. The information is 
received and processed at the second transceiver. The information is differentially encoded using phase 
shift keying. In addition, after a pre-selected time interval, the first transceiver may transmit again. 
During the preselected time interval, the second transceiver may exchange information with another 
transceiver in a time duplex fashion. The processing of the signal at the second transceiver may Include 
estimating the phase differential of the transmitted signal and pre-distorting the transmitted signal. A 
transceiver includes an encoder for encoding information, a wideband frequency division multiplexer 
for multiplexing the information onto wideband frequency voice channels, and a local oscillator for 
upconverting the multiplexed information. The apparatus nnay indude a processor for applying a 
Fourier transform to the multiplexed information to bring the information into the time domain, for 
transmission. 
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FIELD OF THE INVENTION 



This invention relates to voice and data transmission in wireless communications, and particularly between 
fixed and portable transmitters and receivers. 
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CLAIM TO COPYRIGHT 

A portion of the disclosure of this patent document contains material which is subject to copyright protec- 
tion. The copyright owner has no objection to the facsimile reproduction by anyone of the patent document, 
10 as it appears in the Patent Office patent file or records, but otherwise reserves all copyright rights whatsoever! 

BACKGROUND AND SUMMARY OF THE INVENTION 

This patent document presents a new multiple access technique for Personal Communication Networks 
15 (PCN). Personal communication networks are networks that allow individuals and equipment to exchange in- 
formation with each other anywhere at anytime through voice, data or video. PCN typically include a number 
of transceivers, each capable of transmitting and receiving information (voice, dataorvldeo) in theform of elec- 
tromagnetic signals. The transceivers may be fixed or portable, and may be identical or one or more of them 
may be more complex. 

20 The system must allow the transceivers to access each other to enable the exchange of information. When 

there are a number of transceivers, multiple access, that is, access by more than one transceiver to another 

transceiver, must be allowed. 

One of the constraints of designing a PCN is that a transceiver, or portable radio unit, must be small in 

size. The smaller the unit, the better for portability. The small size of the units means only small and light-weight 
25 power sources can be used. If the portable is to be used for any length of time. It must therefore consume 

minimal power. 

Also, to allow use of the radio frequency spectrum without obtaining a license in North America, the system 
must use a spread spectrum and satisfy federal regulations. In part, these regulations impose limits on the 
power and the firequency spread of the signals exchanged between the transceivers. An object of an aspect 

30 of this invention is to satisfy those requirements. 

Also, transceivers talk to each other over a fixed bandwidth. Because of the limited availability of the RF 
spectrum, the system fnust be bandwidth efficient yet at the same time maintain high quality exchange of in- 
formation at all times in one of the most hostile channels known in communication. The new multiple access 
technique proposed here addresses all these issues. 

35 The new access technique has a low Bit Error Probability (BER) as well as a low probability of dropped 

and blocked calls. This is due to the fact that the access technique is robust against multipath, Doppler shifts. 
Impulse noise and narrowband interference. It has a low cochannel interference and little or no intersymboi 
interference. 

The new access technique can offer up to 38 times the capacity of analog FM. It includes in one aspect 
40 wideband orthogonal frequency division multiplexing of the information to be exchanged, and may include slow 
Frequency Hopping (FN). The technique Is implemented using Digital Signal Processors (DSP) replacing con- 
ventional analog devices. The system operates with relatively small cells. In other aspects, dynamic channel 
allocation and voice activation may be used to improve the capacity of the system. 
Advantages of the present invention include: 
45 1 . It can be used indoors as well as outdoors using the same transceivers. If data Is to be exchanged, as 

opposed to voice, the transceiver preferably contains an estimator to allow pre-distortion and post- 
distortion of the transmitted signal. 

2. The system, as compared with prior art systems omits the dock or carrier recovery, automatic gain con- 
trol or passband limiter, power amplifier, an equalizer or an interleaver-deinterleaver, and therefore has 

50 low complexity. 

3. The system offers good speech quality, as well as low probabilities of dropped and blocked calls. It is 
robust against Doppler and multipath shifts, It is also robust against both impulse noise and narrowband 
interference. 

4. The system is flexible, such that at the expense of increased complexity of the DSP receiver it can be 
55 applied over noncontiguous bands. This is accomplished by dividing a 100 MHz (in one of the exemplary 

embodiments described here) band into several subbands each accommodating an integer number of 
voice channels. 

5. The system offers low frame delay (less than 26.2 ms in the exemplary cellular embodiment described 
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here). The transceiver requires low average transmitted power (of the order of 20^W in the exemplary cel- 
lular embodiment described here) which means power saving as well as enhanced biological safety. 
6. The system offers up to a 38 fold increase in capacity over the.North American Advanced Mobile Phone 
System (AMPS) which uses analog frequency modulation. 
5 Operation of the system in accordance with the techniques described in this disclosure may permit conrv 

pliance with technical requirements for spread spectrum systems. 

There is therefore disclosed in one aspectof the invention a method for allowing a number of wireless trans- 
ceivers to exchange information (data, voice or video) with each other. In the method, a first frame of infor- 
mation is multiplexed over a number of frequency bands at a first transceiver, and the information transmitted 
10 to a second transceiver. In a cellular implementation, the second transceiver may be a base station with ca- 
pacity to exchange information with several other transceivers. The information is received and processed at 
the second transceiver. Thefirequency bands are selected to occupy a wideband and are preferably contiguous, 
with the information being differentially encoded using phase shift keying. 

A signal may then be sent from the second transceiver to the first transceiver and de-processed at the 
15 first transceiver. In addition, after a preselected time interval, the first transceiver transmits again. During the 
preselected time interval, the second transceiver may exchange information with another transceiver in a time 
duplex fashion. 

The processing of the signal at the second transceiver may include estimating the phase differential of 
the transmitted signal and pre-distorting the transmitted signal. 

20 The time intervals used by the transceivers nnay be assigned so that a plurality of time intervals are made 

available to the first transceiver for each time interval made available to the second transceiver while the first 
transceiver is transmitting, and for a plurality of time intervals to be made available to the second transceiver 
for each time interval made available to the first transceiver otherwise. Frequencies may also be borrowed 
by one base station from an adjacent base station. Thus if one base station has available a first set of frequen- 

25 cies, and another a second set of distinct frequencies, then a portion of the frequencies in the first set may 
be temporarily re-assigned to the second base station. 

In an implementation of the invention for a local area network, each transceiver may be made identical 
except for its address. 

Apparatus for carrying out the method of the invention is also described here. The basic apparatus is a 
30 transceiver which will include an encoder for encoding information, a wideband frequency division multiplexer 
for multiplexing the information onto wideband frequency voice chiannels, and a local oscillator for upcpnverting 
the multiplexed information. The apparatus may include a processor for applying a Fourier transform to the 
multiplexed information to bring the information into the time domain for transmission. 

35 BRIEF DESCRIPTiON OF THE DRAWINGS 

There will now be described a preferred embodiment of the invention, with reference to the drawings, by 
way of illustration, in which like numerals denote like elements and in which: 

Figures 1a and 1b are schematics of a prjor art receiver and transmitter respectively; 
40 Figure 2 is a schematic showing the use of the available frequencies according to one aspect of the in- 

vention for use with cellular applications; 

Figure 3a is a schematic showing an idealized pulse for transmission over a cellular system; 
Figure 3b is a schematic showing a modified viersion of the pulse shown in Figure 3a; 
Figure 3c is a schematic showing a further modified version of the pulse shown in Figure 3a; 
45 Figure 4 is a scheniatic showing an exemplary protocol for cellular communication; 

Figure 5a is a block diagram showing the structure and function of an embodiment of the transmitter of a 
cellular portable in accordance with the invention; 

Figure 5b is a block diagram showing the structure and function of an embodiment of the transmitter and 
receiver of a cellular base station in accordance with the invention; 
50 Figure 5c is a block diagram showing the structure and function of an embodiment of the receiver of a cel- 

lular portable in accordance with the invention; 

Figure 8a is a flow diagram showing the function of the processor in either of Figures 5a or 5b; 
Figure 6b is a schematic showing the function of the deprocessor in either of Figures 5b or 5c; 
Figure 6c is a schematic further illustrating the operation of the processor and deprocessor shown in Fig- 
55 ures 6a and 6b; 

Figure 7a is a schematic showing the structure and function of the channel estimator in Figure 5b; 
Figure 7b is a flow chart showing the operation of the channel estimator of Figures 5b and 7a; 
Figures 8a, 8b and 8c are respectively schematics of 126, 63 and 7 cell reuse patterns; 
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Figures 9a and 9b are schematics showing transmit protocols according to one aspect of the invention; 
Figure 10 is a schematic showing the use of the available frequencies according to another aspect of the 
invention for use with local area network applications; 

Figure 11a is a schematic showing an idealized pulse for transmission over a local network system; 

Figure lib is a schematic showing a modified version of the pulse shown in Figure 11a; 

Figure 11c is a schematic showing a further modified version of the pulse shown in Figure 11a; 

Figure 12 is a schematic showing a preferred protocol for local area network communication; 

Figure 13a is a block diagram showing the structure and function of an embodiment of the fransmitter of 

a local area network transceiver according to the invention; 

Figure 1 3b is a block diagram showing the structure and function of an embodiment of a further local area 
network transceiver according to the invention; 

Figure 13c is a block diagram showing the structure and function of an embodiment of the receiver of a 
local area network transceiver according to the invention; 

Figure 14a is a flow diagram showing the function of the processor in either of Figures 13a or 13b; 
Figure 14b is a schematic showing the function of the deprocessor in either of Figures 13b or 13c; 
Figure 14cls a schematic further illustrating the operation of the processor and deprocessor shown In Fig- 
ures 14a and 14b; and 

Figure 15 is a schematic showing the structure and function of the channel estimator in Figure 13b. 
DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
Introduction 



The benefits of the invention can be readily appreciated with reference to Figure 1, which shows a prior 
art transmitter/receiver configuration for a portable unit. The transmitter includes a vocoder 110, an interleaver 
112. a modulator 114, a filter 116, local oscillator 118. power amplifier (PA) 120 and antenna 122. The receiver 
includes an LNA 124, a local oscillator 126. a filter 128, automatic gain control (AGC) 130 with an associated 
passband hardllmiter not separately shown, carrier recovery 132, sampler 134, clock recovery 136, adaptive 
(or fixed) equalizer 138, demodulator 140. defnterleaver 142 and decoder 144. With implementation of the 
present invention, several of the blocks shown in Figure 1 are not required. These are the interleaver 112, de- 
interleaver 142, power amplifier 120, automatic gain control 130 with passband hard-limiter, clock recovery 
136 and carrier recovery 1 32. and the equalizer 138. It will now be explained how the proposed system obtains 
the omission of these blocks without impairing the quality and capacity of the system. 

In this disclosure there will be described two systems as examples of the implementation of the Invention. 
The system described first here will apply to a cellular system with a number of portable transceivers and base 
stations (BS). Then will be described a local area network implementation. A local area network will typically 
be a System of equal transceivers. The Invention may also be implemented with combinations of cellular and 
local area networks, or to a system with a number of equal transceivers and a master or controlling transceiver. 
"Equal" as used here means that the transceivers have more or less the same processing equipment and proc- 
essing capabilities. The system described here Is primarily for the exchange of voice Information. 

Link set-up and termination protocols between transceivers, and the equipment required to implement 
them, are well understood in the art as well as the basic structure of radio transceivers that may be used to 
implement the invention. Hence these elements are not described here. What Is described here are the novel 
operational, functional and structural elements that constitute the invention. 

Cellular Implementation of Wideband Modulation 

The present invention proposes In one embodiment a wideband modulation scheme for exchange of In- 
formation between transceivers such as portables and base stations. 

Wideband in this patent document is described in the context of Wideband-Orthogonal Frequency Domain 
Modulation (W-OFDM or wideband OFDM). In OFDM, the entire available bandwidth B Is divided into a number 
of points K, where adjacent points are separated by a frequency band Af, that is B = KAf. The K points are 
grouped into a frame of Ki points and two tail slots of K2 points each, so that K = Ki + 2K2. The frame carries 
the information Intended for transmission under theform of multilevel differential phase shift keying (MDPSK) 
symbols or differential quadrature amplitude modulated (DQAM) symbols. Thus each point in the frame cor- 
responds to one information symbol. The two tail slots act as guard bands to ensure that the out-of-band signal 
is below a certain power level. For example, when a pulse P(f) is selected for pulse shaping and the out-of- 
band signal has to be ydB or less relative to the in-band signal. K2 is selected such that 
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20.logio I P(f)/P(6) ( ^ y for f ^ KjAf: 
When the pulse is a raised-cosine pulse with a roll-off p and when the number of levels each symbol can take 
Is M, the bit rate Is equal to Kilog2M/(6t + (1 + p)/Af) where (1 + p)/Af is the duration of the frame and.6t is the 
guard time required to take into account the delay of arrival and the delay spread due to multlpath. In this case, 
5 the bandwidth efficiency, which is defined as the ratio between the bit rate and the bandwidth, is equal to: 

log2M/((1 + p + 6tAf)(1 + 2K2/K1)) 
In wideband-OFDM, both K and Af are selected sufficiently large to achieve a high throughput as well as 
to reduce the effects on the BER of the clock error, the Doppier shift and the frequency offset between the 
LO in the transmitter and the one in the receiver. To show what is meant by "K and Af aire selected sufficiently 
- 16 large", consider the effect of increasing K and Af on (1) the clock error, (2) the Doppier shift and (3) the fre- 
quency offset between the LO in the transmitter and the LO in the receiver. 

(1) When a dock error at a transceiver of value x occurs in the time domain, it causes a shift in the phase 
difference between adjacent symbols In the frequency domain of value 2nMx, When x is equal to xT where 
T is duration of one time domain sample and % is any real value, the shift is equal to 27cAfxT. Hence, x causes 

15 a shift in the phase difference between adjacent symbols of value 2nxlKy since T is equal to 1/(KiAf). By 

doubling the number of symbols from Ki to 2Ki the shift in the phase difference is reduced by half from 
2n%IK^ to Tix/Ki- Thus, the effect of the clock error on the BER is reduced by increasing K. 

(2) When there is relative motion between the transmitting transceiver and the receiving transceiver, a Dop- 
pier shift occurs with a maximunri absolute value I V/x| where V is the relative velocity between the two 

20 transceivers and X is the wavelength of the travelling wave corresponding to the carrier frequency f^ (i.e. 

fc is the frequency corresponding to the middle point in the frame). Such a Doppier shift causes a sampling 
error in the frequency domain of the same amount, or equivalently. it causes a sampling error of V/(XAf) 
relative to one symbol sample. Thus, the effect of the Doppter shift on the BER is reduced by increasing 
Af. 

25 (3) When a frequency offset between the LO in the transmitter and the one In the receiver occurs with a 

value fo, it causes a sampling error in the frequency domain of the same amount, or equivalently, it causes 
a sampling error of fo/Af relative to one symbol sample. Thus, the effect on the BER of the frequency offset 
between the LO in transmitter and the one in the receiver is reduced by increasing Af. 
In summary, OFDM with a K and a Af large enough to be able to achieve a specific throughput and large 
30 enough to be able to avoid using either a dock or a carrier recovery device without substantially affecting the 
BER is referred to here as Wkjeband-OFDM. As an example, let us assume that MDPSK is used in an OFDM 
system with the number M of levels, with a carrier frequency f^, with a raised cosine pulse of roll-off p, with 
the LO at the receiver having a frequency offset fo relative to the LO at the transmitter (so that the frequency 
offset between the carrier frequencies in the first and second transceivers of the multiplexed Information is 
35 fo), with a given maximum expected clock error x = xT at the receiving transceiver, where T is the duration of 
one time domain sample, and with a maximum expected relative velocity V between the transceivers. Thus, 
in order to ensure that the out-of-band signal is ydB or less relative to the in-band signal and to be able to 
avoid using either a clock or a carrier recovery device without substantially affecting the BER we have to: 
1.. Find the acceptable sampling error Af, relative to one symbol sample, which does not substantially affect 
40 the BER. This can be done using the following rules: 

When 0.2 ^ p ^ 0.3, Af = 7.50% 
WhenO.3 ^P ^0.4, Af = 10.0% 
When 0.4 ^ p ^ 0.5, Af = 12.5% 
When 0.5 ^ p ^ 0.6, Af = 15.0% 

45 2. Find Af such that: 

V/(XAf) + fo/Af ^ Af ' 

3. Find K2 such that 

20.logio I P(f)/P(0) I ^ y for f ^ K2Af 

4. Find Ki such that 

50 2nx/Ki < nM 

In this case, we refer to OFDM as Wideband-OFDM. Element 4 is a necessary condition for wideband OFDM, 
and given a sampling error, the sampling error may be corrected with the methods described in this patent 
document 

To implement wideband modulation, Orthogonal Frequency Division Multiplexing (OFDM)) is preferred in 
55 which the inforrhation, for example encoded speech, is multiplexed over a number of contiguous frequency 
bands. Wideband OFDM forces the channel to be frequency selective and causes two types of linear distortion: 
amplitude distortion and phase distortion. To reduce the effect of amplitude distortion the modulation is pre- 
ferably phase modulation, while the effect of phase distortion is reduced by employing differential phase mod- 
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ulation. Hence the modulation may be referred to as Differentia! OFDM (DOFDM). Unlike in other proposed 
schemes, neither pilot tones nor diversity are required in DOFDM. Possibly, quadrature amplitude modulation 
might be used, but amplitude modulation makes it difficult to equalize the distorting effects of the channel on 
the signal. 

5 To implement wideband modulation in a cellular system with a plurality of portables and one or more base 
stations, a 100MHz band is divided into 4096 points, as shown in Figure 2, plus two tail slots of 195.3 KHz 
each. The 4096 points represent N voice channels (vc). Adjacent points are separated by 24.414 KHz and each 
point represents a Differential eight Phase Shift Keying (D8PSK) Symbol eJ«n), where ^(n) = C(n-1) + 0(n) + 
X(n). 0(n) takes one of the eight values {0, 27i/8, 47c/8 147r/8} with equal probability for n = 1, 2, 4096 

io and 0(0) takes an arbitrary value. x(n) also takes an arbitrary value. x(n) may be used as a security key and 
will be known only to the transmitter and receiver. Information in the form of output bits from a vocoder are 
mapped onto 0(n). Vocoders are well known in the art and do not need to be described in detail here. The focus 
here is to transmit the bits with an acceptable Bit Error Rate, I.e. with a BER ^10-2 for voice and ^10-» for 
data. 

t5 Ideally, 40.96 ^ls (=1/24.414 KHz) is the minimum duration required for one frame to be transmitted without 

firequency domain rntersymbol interference. This can be achieved using a Raised Cosine (RC) pulse with zero 
roll-off, as shown in Figure 3a. Figure 3a illustrates a rectangular (time domain) window corresponding to the 
RC (frequency domain) pulse. Such a pulse, however, requires an infinite frequency band. To alleviate such a 
requirement, an RC pulse with a 20% roll-off (i.e. p = 0.2) may be used as shown in Figure 3b. The frame dur- 

20 ation has increased by 20% to 49,1 52 us. The two tail slots of 195.3 KHz each (i.e. 8 points each) ensure that 
the signal outside the entire band of 100.39 MHz is below -50 dB. To allow the frame to spread over the time 
as a consequence of the multipath nature of the channel, an excess frame duration of 2.848 ^s is provided as 
shown in Figure 3c, making the frame duration 52 jis in total. 

Since the frame duration is 52 jas, the frame rate is 252 frames per 13.1 04 ms or equivalently, 126 full du- 

25 plex frames may be transmitted/received every 13.104 ms. The reason for pre-selecting an interval of 13.104 
ms is to ensure a transmission delay to allow one transceiver to communicate with other transceivers at the 
same time, but must not be so long that the delay becomes unacceptable to the user. Delays longer than about 
40 ms are too great for voice, and it is preferable to be lower. For data, the delay may be longer and still be 
acceptable. 

30 In the exemplary embodiment described here, three bit rates are considered for the vocoder: 18.77 Kbps, 

9.16. Kbps and 6.18 Kbps. Table I displays the structure of a vc slot and the number N of vc for each vocoder 
rate. The control symbols in each vc slot are required for handoff and power control. Figure 2 shows that N vc 
can be transmitted simultaneously. This Is known as Frequency Division Multiple Access. Figure 3c shows that 
126 full duplex frames can be transmitted every 13.1 04 ms in a Time Division Multiple Access fashion (TDMA). 

35 The total number of Full Duplex voice channels (FDvc) is therefore 126xN and is shown in Table 1. 

To ensure that the channel is slowly fading, a Time Division Duplex protocol for exchange of information 
between the portable and the base station is proposed as illustrated in Figure 4. The protocol is as follows: 

1 . The portable transmits a frame 410 over one vc slot. See the discussion in relation to Figure 5a below. 

2. The Base Station (BS) receives the frame 410 from the portable and processes (analyzes) it as shown 
40 and discussed in relation to Figure 5b below. 

3. Based on the received signal, the BS predistorts a frame 420 and transmits it to the portable over the 
same vc slot, 520 ^s or some other suitable time interval later in which the channel does not change sub- 
stantially. The time interval will depend on factors such as the frequency, speed of the fransceiver and 
other environmental factors. 

45 4. The portable receives the frame from the BS. See the discussion in relation to Figure 5c below. 

5. Steps 1 through 4 are repeated, as for example by the transmission of the next frame 430, every 13.104 
ms until the call is terminated. 

During 520 ^is, a portable travelling outdoor at 1 00 km/hr moves 1 .44 cm, which leaves the outdoor channel 
largely unchanged. Indoors, a portable moving at 2 m/s moves 0.1 cm again leaving the channel unchanged. 
50 Assuming that the channel is reciprocal and stationary over 520 us, a predistorted signal, transmitted by the 
BS, should reach the portable undistorted. 

From Figure 4, one can see that the portable transmits/receives one FDvc every 13.104 ms, while the BS 
can transmit/receive up to 21 frames or equivalently up to 21 xN FDvc every 1 3.104 ms. The frames 440 labelled 
frame 2.. .frame 21 are frames that may be transmitted to other portables. This implies that while one BS is 
55 processing its data over 520 |as, six other BS can communicate to their corresponding portables in a Time Di- 
vision Multiple Access (TDMA) fashion using the same frequency bands. Also, during the 13.104 ms, or such 
other preselected time interval that is suitable, the BS may communicate with one or more other portables. 

When a portable is stationary during a call, it is possible with high probability to have the transmitted signal 
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centered with several deep (frequency domain) nulls, hence, causing speech degradation. Also, narrowband 
interference over the vc slot can deteriorate the speech. In order to avoid both situations, the signal is prefer- 
ably frequency hopped into a new vc slot within the same (frequency domain) frame. This frequency hopping 
is ordered by the BS which is constantly monitoring the channel frequency response. Monitoring techniques. 

5 as well as frequency hopping, are known in the art, and not described here further. When an unacceptable 
speech degradation is first noticed by t he BS a probation period is initiated and maintained for at least 1 0 cycles 
(i.e. 10x13.104 ms) unless speech degradation has ceased. In other words, the probation period is terminated 
if speech degradation has ceased. Frequency hopping is then ordered at the end of the probation period. The 
period of 10 cycles is long enough to Indicate the portable statlonarity and is short enough to allow speech 

10 interpolation between unacceptable speech frames, hence maintaining good speech quality. As known In the 
art, the BS ensures that no collisions take place between hopping portables. 

Digital Signai Processing 

15 The transmitter/receiver block diagrams corresponding to the protocol In Figure 4 are shown in Figures 

5a, 5b and 5c. Figure 5a corresponds to step 1 in the protocol described above. Speech is provided to a vocoder 
510 where the speech is digitized and coded to create bits of information. The bits are provided to the modulator 
51 2 which turns them intoDSPSKsymbols, with three bits persymbol.The D8PSK symbols are then processed 
in the processor 514 which is described in more detail in Figure 6a. The output from the processor is then fil- 
20 tered in low pass filter 516, upconverted to RF frequencies using local oscillator 51 8 and transmitted by antenna 
520. Figure 5b corresponds to steps 2 and 3. 

In Figure 5b, the received signal at the base station is filtered In a bandpass filter 522, and down converted 
by mixing with the output of a local oscillator 524. The average power of the downcoverted signal is monitored 
by a power controller 525 that adjusts the average power to the specifications required by the sampler 526. 
25 The adjusted downconverted bits are then sampled in sampler 526 to produce bits of Information. The bits are 
then processed in the deprocessor 528, described in more detail in Figure 6b. An estimate of the phase dif- 
ferential of the received signal Is taken in the channel estimator 530, as described in more detail in relation to 
Figure 7a and 7b below, and the estimated phase differential Is supplied to a decoder-demodulator 532 to cor- 
rect the received bits. The estimated phase differential is also supplied to a pre-distorter 534 in the transmitter. 
30 At the transmitter in the Base Station, the same blocks are Incorporated as in the portable transmitter except 
that a pre-distorter is used to altier the phase of the D8PSK symbols to make the channel appear Gaussian 
(ideal) as opposed to a fading channel. The pre-distorter 534 receives a signal corresponding to the estimated 
phase differential of the channel. On the (believed reasonable) assumption that the channel is reciprocal, the 
signal being transmitted is predistorted with the estimated phase differential so that the received signal at the 
35 portable with which the BS is communicating will be corrected for any phase distortion over the channel. The 
advantage of rendering the channel Gaussian is a large saving in the power required to achieve an acceptable 
BER. The initial power control 525 also sends a signal to the pre-distorter 534 to adjust the transmitted power 
to an appropriate signal level for the sampler 526 in the portable's receiver depending on the average power 
of the received signal. Thus if the average power is too low, the transmitted power is increased and if the aver- 
se age power is too high, the transmitted power is decreased. The power controller 525 may also be used in fre- 
quency hopping to monitor the average power of the received frequency and determine when frequency hop- 
ping need take place. 

Figure 5c corresponds to step 4, and shows the receiver of the portable, which is the same as the receiver 
in the BS except it does not include an estimator or a power controller. These are not required in the portable 

45 on the assumption that the BS will carry out the phase estimation and the power control. However, if desired, 
the portable may include these functions. 

Figures 6a, 6b and 6c illustrate the function and. structure of the processor and the deprocessor respec- 
tively in the transmitter and receiver. Software for modelling the function of the processor in a general purpose 
computer is attached hereto as Appendix A, and for modelling the function of the deprocessor as Appendix B. 

50 Figure 6a shows that the processor is a DSP implementation of an RC pulse shaping filter with a 20% roll- 

off, followed by an inverse Fourier transform. The processor first inverse Fourier transforms the 4096 D8PSK 
modulated synibols output from the modulator. The transformed symbols are then triplicated as a group so 
that the total number of samples Is tripled, with three consecutive groups each consisting of the 4096 trans- 
formed symbols. The triplication of the signal is illustrated in Figure 6c, where the symbols are shown as first 

55 delayed and added together. Next, as shown in Figures 6a and 6c, the three groups are windowed by a Raised 
Cosine window with a roll-off of 0.2 centered in the middle of the three groups. In other words, the processor 
takes D8PSK symbols in, pulse shapes them and inverse Fourier transforms them. On the other hand, the 
deprocessor undoes what the processor did, i.e. it removes the pulse shaping, then Fourier transforms the 
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received signal to obtain llie original D8PSK symbols. The first two blocks in Figure 6b are similar to the second 
two blocks in Figure 6a except for two differences. The two differences are as follows. In the first block of the 
deprocessor, the repeated groups of symbols are partially overlapped as shown in Figure 6c on the right hand 
side. In the second block, a rectangular window is used instead of the Raised Cosine. In the preferred imple- 
mentation, the blocks are repeated three times but other numbers of repetition may be used. 

Figures 6a, 6b and 6c show that the DSP blocks used In the processor are Identical to the ones used In 
the deprocessor, except for a small change in the two transforms and a small change in the shapes of the two 
windows. Thus the same hardware can be used by both the processor and the deprocessor. 

Figure 7a shows a block diagram of an example of a preferred channel estimator, and Figure 7b is a flow 
chart showing the operation of the phase estimator. Each of the steps is carried out in a computing means 
that may be a special purpose computer or a general purpose computer programmed to carry out the digital 
signal processing described here, as for example with the software that is attached hereto as appendix C. 
Other methods of estimating the channel may be used that obtain an estimate of the channel group delay or 
phase differential of the transmitted symbols. However, a preferred implementation Is described here. 

The f irst block in Figure 7a estimates the envelope A(n) for n = 1 , 4096 of the (frequency domain) sam- 
ples transmitted over the fading channel as output from the deprocessor, which may be filtered in accordance 
with known techniques before or after estimation of the envelope. The estimate A*(n) is the square-root of the 
sum of the squares of the quadrature (Q) and inphase (I) samples output from the deprocessor. The second 

block performs the operation: Aln(A'(t)) = (A^n) - A(n-1))/A(n), for n = 2 4098, where A'(n) is the estimate 

of A(n). The third block performs a Hilbert transform operation H[Aln(A'{t))] on the result of the second block. 

HlAln(A'(t))] is an estimate of | Aco(n) I for n = 2 4096, where A(D(n) is the phase differential of the transmitted 

signal (© is the phase of the signal). The Hilbert transform is preferably carried out by taking the discrete fast 
Fourier transform of the data record, multiplying the positive frequency spectrum of the transform by-i (square 
root -1 ), and the negative frequency spectrum of the transform by i, and taking the inverse discrete fast Fourier 
transform. The result is a set of symbols representing an estimate of the phase differential of the received sig- 
nal, as determined from its sampled amplitude envelope. 

Instead of a Hilbert transform, a different estimation may be made to estimate the phase differential. In 
this case, f irstly, after the evelope and electromagnetic signal has been sampled, a series of data frames of 
a number of consecutive amplitude samples (A{t)) of the electromagnetic signal are constructed. These data 
frames are then segmented into segments [UM. where the amplitude of the electromagnetic signal is at least 
a predetermined number of dB less than its running mean, for example, 10dB. The following calculation is then 
applied to these segments of the amplitude samples: 

A(D(t) — — — - 

where t' = t - t^fn. t^in is the time in [ti, tj when A(t) reaches its minimum, t is the time from the beginning of 
the segment, and t© is the time from the instant the amplitude of the electromagnetic signal reaches its minimum 
during the segment until the amplitude reaches double its minimum during the segment. In other words, the 
phase differential may be calculated from 

Ac£)(t) « - to/(to2 + r2). 

The polarity of A(D(n) is extracted using the last block shown in Figure 7a. The estimate so calculated does 
not provide the sign of the differential. This may be determined by known techniques, for example by adding 
the phase differential to and subtracting the phase differential from the received phase (tan-^ (Q/l)) and taking 
the sign to be positive if the addition results In the smaller Euclidean distance to the expected value and neg- 
ative if the subtraction results in the smaller Euclidean distance to the expected value. 

Equivalently, for each sample n, the ideal phase closest to (o(n)+AcD(n) is determined and labelled Mn), 
and the ideal phase closest to co(n)-Aco(n) is determined and labelled cD_(n). The two sums P = 2:|<»4<n) - 
{a)(n)+AcD(n)} I and N = 5:|o>__(n) - {(D(n)-Aco(n)} | are calculated. If P < N, then o>(n)+A(D(n) is used to correct 
the signal, and if not then G)(n)-A(D(n) is used to correct the signal. 

For simplicity of the estimator, the determination of the sign need only be carried out for phase differentials 
greater than a predetermined threshold. This will be in the vicinity of a fade and may be accomplished by seg- 
menting the data record into a segment In which the phase differential is larger than a selected threshold and 
setting the remainder of the data record to zero. This computation may be carried out with a simple discrim- 
ination circuit or equivalent computing means in the estimator. 

The bias 6cd of the channel group delay Is estimated by averaging A(o'(n) over n for n = 1, 4096 where 
AQ)'(n) is the measured value of Aco(n). The estimates A'(n) and Aa)'(n) are used directiy In the predistortion 
filter in Figure 5b. while the estimates A(o{n) and 5© of the unbiased channel group delay and of the bias of 
the channel group delay respectively are used in the demodulator. 
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The complexity of the processor-deprocessor-channel estimator is displayed in Table II. Complexity is 
measured in Mega Instructions Per Second (MIPS) where one instruction Is defined as one complex addition, 
one complex multiplication and a storage of one complex number. It does not include overhead. 

The complexity of the processor-deprocessor-channel estimator in the BS is computed firom the complex- 
5 Ity of the Inverse Fast Fourier Transform (IFFT)/Fast Fourier Transform (FFT)/Hilbert Transform. The complex- 
ity is 4096 X 12 X 4 X 21 /13.104ms for the BS. For the portable, it is computed from the complexity of the FFT- 
/IFFT per vc: (32 x 5 + 64 + 128 + 256 + 512 + 1024 + 2048 + 4096)2/1 3.1 04ms for the portable with a 6.18 
Kbps vocoder. Such a complexity assumes that the A/D converter operates at 100 MHz with 12 bit precision. 
As seen in Table II, the portable has smaller complexity due to the fact that the portable transmits/receives 
10 one vc in 13.104 ms and the BS transmits/receives up to 21xN vc in 13.104 ms. 

Reducing Analog Complexity 

Comparing Figure 1 (prior art) and Figure 5. it will be seen that several conventional blocks are not used 
IS in the present invention, namely the interleaver-deinterleaver. the Power Amplifier (PA), both the clocl< and 
the carrier recovery, both the AGC with its associated Passband hard limiter, as well as the equalizer. 

From the BS point of view, the interleaver-deinterleaver is not required since the signal is predistorted 
before transmission forcing the received samples to be independent. From the portable point of view, the in- 
terleaver-deinterleaver is not required as a separate entity from the vocoder due to the fact that the channel 
20 is highly frequency selective, hence the interleaving/deinterleaving can be applied implicitly in the vocoder over 
one vc, without a need for a separate time domain tnterleaver/deinterleaver. This eliminates excess speech 
delays associated with interleaving/deinterleaving between frames. 

The PA is not required since the cells can have, as shown later, a radius of up to at least 250 m outdoors 
and 30 m indoors, if the transmitted power is up to 6 dBm. Such a power can be generated by the Local Os- 
25 clllator (LO) without a need for a PA. It is important to avoid using a PA since DOFDM generates a time domain 
signal with non constant envelope. A power efficient class C PA cannot be used without distortirig the signal. 
A dass A PA can be used at the expense of power efficiency. 

A dock recovery device is not required since a sampling error in the time domain is equivalent to a phase 
shift in the frequency domain. The phase shift is a linear function of frequency. It contributes to the bias in the 
30 . channel group delay. Such a bias can be easily estimated and removed as mentioned previously by averaging . 
o'(n) over n in the frequency domain. Such an estimate is accurate as long as the sampling error is. less than 
0.2 \iS or equivalently less than 20 samples (since in this case, the corresponding phase shift is less than n), 
and as long as the number of points in one vc is large enough as it is here. 

Acarrier recovery device is not required since a carrier offset in the time domain is equivalent to a sampling 
35 error in the frequency domain. For the chosen RC pulse, a sampling errorof up to 10% of the duration of one 
pulse is acceptable. This implies that a frequency offset of up to 2.414 KHz is acceptable regardless whether 
it is due to carrier offset as low as 1 part in a million, i.e. as low as 1 KHz per 1 GHz. When a carrier frequency 
higher than 2.414 GHz is required, one can decrease in Figure 2 the numbei* of points per 100 MHz or one 
can use an RC pulse with a rolloff larger than 20%. 
40 Neither an AGC nor a Passband hard-limiter are required since the level of the received power may be 

controlled constantly. This is achieved as follows: The portable transmits a frame. The BS receives the frame 
and predistorts a frame intended for transmission accordingly, assuming that the channel is reciprocal and sta- 
tionary over 520 ^is. This indudes controlling the transmitted power according to the received power. The BS 
transmitsthe predistorted frame and simultaneously orders the portable to control its power. The order is con- 
45 veyed using the control symbol in the vc slot (See table I). The degree of power control may be determined 
using the power controller 525, and the instruction for the inclusion of a power control symbol in the vc may 
be sent from the power controller 525 to the pre-distorter 534. 

One advantage of wideband modulation over narrowband modulation is that the wideband signal does not 
experience short term fading the same way the narrowband one does. The wideband signal is mainly affected 
50 by shadowing and other long term effects which vary slowly and are easily monitored from one frame to the 
other as long as the same vc slot Is used by the portable to transmit and receive (i.e. as long as TDD is em- 
ployed). 

Finally, conventional equalization, whether it is linear or nonlinear, is not required simply because there is 
little or no ISI. Also, frorn the portable point of view, each received vc is predistorted by the BS. Hence, the 
55 channel can be modeled approximately as an ideal memoryless Additive White Noise Gaussian (AWGN) chan- 
nel, assuming channel reciprocity and stationarity over 520 ^s. From the BS point of view. since the received 
signal is not predistorted by the portable prior transmission, the channel estimator is used to reduce the effect 
of the channel group delay. 
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Smaller cells 

As mentioned previously, the LO generates a 6 dBm average power, hence the signal power transmitted 
by the BS over one vc slot is (6 dBm - lOlogioN dB) while the signal power transmitted by the portable over 
5 one vc slot is 0 dBm. Also, since the noise power over a 100 MHz band is -94 dBm, ft is (-94 dBm - lOlogioN 
dB) over one vc. Atypical noise figure at the receiver is 7 dB. The penalty for not using a matched filter in the 
receiver is 1 dB. Combining together the above figures provides the portable with an (92 dB - path loss in dB) 
received signal to noise ratio (SNR), while it provides the BS wfth an (86 dB + lOlogipN dB - path loss in dB) 
received SNR. 

10 For a path loss of 75 dB. the radius of the urban ceil can be 250 m while it can be 30 m for the indoor cell. 

Such a path loss provides the portable wfth a 17 dB received SNR, while it provides the BS with an (11 dB + 
lOlogioN dB) received SNR. From the portable point of view, the channel can be modeled approximately as 
an ideal AWGN channel, hence the 17 dB received SNR results in a 2x10-^ BER. On the other hand, the chan- 
nel can be pessimistically modeled as a Rayleigh fading channel from the BS point of view. The corresponding 

IS BER are displayed in Table III which shows that the achieved BER is ^ 4x10-^. A BER ^ IO-2 is acceptable 
for voice. 



Cell Pattern Reuse 

20 From Table I, the number of Full Duplex voice channels (FDvc) that can be transmitted/receh/ed per frame 

is 136 over 100 MHz, for a 6.18 Kbps vocoder. If the bandwidth is halved to 50 MHz. the number of FDvc per 
frame is reduced to 68. the noise floor is reduced by 3 dB and the number of full duplex frames that a BS can 
transmit/receive is doubled to 42. leaving the frame duration, the number of frames per 13.104 ms and the 
processor/deprocessor complexfty unchanged. 

25 Reducing the available bandwidth directly affects the cell pattern reuse. This can be explained as follows, 

assuming that we are required to offer a minimum of 1 36 FDvc per cell, that the vocoder rate is 6.18 Kbps and 
that the cell radius is fixed at 250 m outdoors and 30 m indoors. For a 100 MHz band, we assign one frame 
per cell and offer 1 36 FDvc per cell. In this case, the cell pattern reuse consists of 126 cells as shown in Figure 
8a which displays a seven layer structure. For a 50 MHz band, we assign two frames per cell and offer 136 

30 FDvc per cell, hence reducing our cell pattern reuse to a 63 cell pattern as shown in Figure 8b which displays 
a five layer structure. If the available bandwidth is as low as 5.86 MHz. we have 8 vc per frame. Hence we 
have to assign 1 8 frames per cell In order to offer the minimum required number of FDvc per cell. This reduces 
the cell pattern reuse to as low as a 7 cell pattern as shown in Figure 8c which displays a two layer structure. 
In Figures 8a, b and c, a shaded area is shown around the center of the pattern, indicating 19, 38 and 126 

35 ftjil duplex frames that the central BS can transmit/receive respectively. Tables IVa, b and c show the number 
of cell layers In each cell pattern reuse, the coverage area in Km2 of the pattern reuse for both the indoor and 
the urban environments, as well as the carrier to interference ratio (CIR) in dB, for the 1 00 MHz. 50 MHz and 
5.96 MHz bands, respectively. In all cases, the CIR is large enough to sustain a toll quality speech. 



40 Transmission/Reception Protocol 

Since the number of FDvc a portable can transmit/receive is one, while the number of FDvc a BS can trans- 
mit/receive is much larger as shown in Table V for each of the three vocoder rates, we have chosen the following 
transmission/reception protocol: 
45 1 . The portable transmits a frame over a vc. 

2. Seven adjacent BS receive the frame from the portable. 

3. One BS transmits to the portable, depending for example on the strength of the received signal by each 
of the BS. 

The control of this protocol may use any of several known techniques. For example, the commonly used 
50 technique is to have the portable monitor the channel and determine which of several base stations ft is closest 
to. It can then order the nearest BS to communicate with it. Another technique is to use a master control which 
receives information about the strength of the signal on the channel used by the portable and controls the 
BS accordingly. Such techniques in themselves are known and do not form part of the invention. 

Such a protocol has several advantages. For instance, the location of the portable can be determined with 
55 high accuracy based on the received vc at the seven adjacent BS. Locating the portable can assist In the BS 
hand-off. A BS hand-off and a portable hand-off do not necessarily occur simultaneously, contrary to other 
prior art systems. In the present invention, when a portable roams from one cell X to an adjacent cell Y, a new 
vc is not required immediately. What is required is a BS hand-off, meaning that BS Y (associated with cell Y) 
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must initiate transmission to the portable over the same vc, while the BS X (associated with cell X) must ter- 
minate its transmission to the portable. 

A BS hand-off occurs without the knowledge of the portable and can occur several times before a portable 
hand-off is required. A portable hand-off is required only when the CIR is below a certain level. In this case. 
5 the Mobile Telephone Switching Office (not shown) calls for a portable hand-off in accordance with known pro- 
cedures. Reducing the portable hand-off rate reduces the probability of dropped calls. This is because a drop- 
ped call occurs either because the portable hand-off is not successful or because there are no available chan- 
nels in cell Y. . 

The present invention allows the use of post-detection diversity at the BS. and the use of dynamic channel 
10 allocation (DCA). 

Dynamic Channel Allocation 

DCA is made possible by each BS having capability to transmit/receive more than the number of FDvc 

IS allocated to Its cell, namely seven times the number of FDvc for a 5.86 MHz band and up to twenty-one times 
the number of FDvc for a 100 MHz as well as a 50 MHz band. The DCA protocol simply consists of borrowing 
as many FDvc as needed from the adjacent cells, up to a certain limit. The limit for the case when we employ 
a 6.18 Kbps vocoder, a 5.86 MHz band and 18 frames per cell is obtained as follows. The cell reuse pattern 
consists of 7 cells. Each cell is preassigned 144 FDvc, Assuming that at peak hours. 75 FDvc are used on the 

20 average and 5 FDvc are reserved at all times, then we are left with 64 idle channels which represent the limit 
on the number of FDvc one can borrow from the cell. 

One should distinguish between the limit on the channels borrowed and the limit on the nonpreassigned 
channels a BS can use. For instance, If a call originates in cell X and the portable roams into an adjacent cell 
Y where no preassigned cells are available. BS Y does not need to borrow immediately a new channel from 

25 an adjacent cell, it can use the original channel as long as the level of CIR is acceptable. If on the other hand, 
a portable wants to initiate a call in cell Y where all preassigned channels are used, BS Y can borrow a channel 
from an adjacent cell up to a limit of 64 channels per cell. 

The main advantage of DCA over Fixed Channel Allocation (FCA) is the increase in traffic handling capa- 
bility. For FCA, a 7 cell pattern each with a preassigned 144 Fdvc can carry a total traffic of 880.81 Eriang at 

30 0.01 Blocking Probability (BP). For DCA. a 7 cell pattern consists of 6 cells each with,80 FDvc that can carry 
a total traffic of 392.17 Eriang, combined with one cell with 528 FDvc that can carry 501 .74 Eriang. The total 
traffic is therefore 893.91 Eriang. This increase appears to be marginal (1.5%). However, if 501.74 Eriang are 
actually offered to one cell in the FCA system (with 14 FDvc/cell), while the six other cells carry 392.17/6 = 
65.36 Eriang per cell, the BP at that busy cell 0.714 while it is negligible at the six other cells. The total blocked 

35 traffic (i.e. lost traffic) in the FCA system is then equal to (6x65.36x0.0 + 1x0.714x501.24) = 358.24 Eriang. 
This represents a 0.4 average BP. If the DCA Is allowed such a loss, its traffic handling capacity would increase 
to 1768.04 Eriang which represents a 100% increase in traffic handling capacity over the FCAsystem, or equiv- 
alently a 160% increase in the number of available FDvc. The DCA system thus represents a marked improve- 
ment over the FCA system. 

40 

Voice Activation 

Voice activation is controlled by the BS according to techniques known in the art. At any instant during a 
conversation between a BS and a portable, there are four possibilities: 
45 1. BS talks while the portable listens. 

2. BS listens while the portable talks. 

3. BS and portable talk simultaneously. 

4. BS and portable listen simultaneously. 

The BS controls the voice activation procedure by allocating in cases 1, 3 and 4 three slots (frames 1.1. 
50 1 .2 and 1 .3) to the BS and one slot the portable (frame 1) every four slots as shown in Figure 9a, Likewise up 
to 21 portables may communicate with the base station in like fashion. 

In case 2, on receiving a signal from the portable, the BS allocates three slots (frames 1.1, 1.2 and 1.3) 
to the portable and one slot (frame 1) to the BS every four slots as shown in Figure 9b. Likewise, up to 21 
other portables may communicate with the base station in like fashion. Consequently, instead of transmitting 
55 two full duplex voice frames over four slots as in Figure 4, voice activation allows us to transmit three full duplex 
voice frames over four slots. Hence, voice activation provides a 50% increase in the number of available-FDvc 
at the expense of increasing DSP complexity. 
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Capacity 

The capacity of Code Division Multiple Access (CDMA) may be defined as the number of half duplex voice 
channels (HDvc) effectively available over a 1.25 MHz band per cell. Based on such a definition. Table IV dis- 
5 plays the capacity of analog FM and of the present system with a 6.1 8 Kbps vocoder. 5.86 MHz band, 1 frame 
per cell and DCA. As shown in Table IV, the capacity of analog FM is 6 HDvc/1 .25MHz/cell while for the present 
system it Is 150 HDvc/1 .25MHz/celL The 6.25 MHz band consists of 5.86 MHz plus two tail slots. When voice 
activation Is used, the capacity of the present system is increased by 1 .5 times to 225 HDvc/1 .25MHz/cell, a 
38 fold increase over analog FM. 

10 

Local Area Networks 

The Invention may also be applied to produce a 48 Mbps wireless LAN, which also satisfies the technical 
requirements for spread spectrum. 

15 For wireless LAN, wideband differential orthogonal frequency division multiplexing is again employed. The 

LAN will Incorporate a plurality of transceivers, ail more or less equal in terms of processing complexity, and 
possibly with identical components except for addresses. 

To implement wideband modulation for a LAN, a 26 MHz band is divided into 1 28 points, as shown in Figure 
10, plus two tail slots of 1.48 MHz each within the 26 MHz band. Adjacent points are separated by 180 KHz 

20 and each point, as with the application described above for a portable-base station, represents a D8PSK sym- 
bol. The transmitter components will be the same as shown in Figure 5b, with suitable modifications as de- 
scribed In the following, and will Include an encoder. The output bits from the encoder are mapped onto the 
D8PSK symbols. 

The frame duration for the symbols is illustrated in Figure 11. A rectangular time domain window corre- 
25 spending to a RC frequency domain pulse has a 5.55 us duration, and includes a 25% roll-off and excess frame 
duration of 0.26 us, making a total 7.2 us duration for the frame. 

For such a wireless local area network (LAN), in which the transceivers are equal, the Time Division Duplex 
protocol is as Illustrated in Figure 12 (assuming there are at least a pair of transceivers): 
1 . A first transceiver transmits a signal (frame 0) over the entire frame. 
30 2. A second transceiver receives the signal from the first transceiver and processes (analyzes) It. 

3. Based on the received signal, the second transceiver predistorts and transmits nine frames (frames 1 
- 9) to the first transceiver immediately. 

Each transceiver has transmitter components similar to those illustrated in Figure 5b. with suitable modi- 
fications to the internal structure to allow the use of the particular frequency band and frame duration em- 
35 ployed. 

The transmitter/receiver functional and sfructural block diagrams are shown in Figures 13a, 13b and 13c 
for the exchange of data. Data is provided to an encoder 810 where the data is digitized and coded to create 
bits of Informatwn. The bits are provided to the modulator 812 which turns them into D8PSK symbols, with 
three bits per symbol. The D8PSK symbols are then processed in the processor 814 which is described in more 

40 detail in Figure 14a. The output from the processor Is then filtered in low pass filter 816, upconverted to RF 
frequencies using local oscillator 818 and transmitted by antenna 820. 

In Figure 13b, the received signal at the base station is filtered in a bandpass filter 822, and down con- 
verted by mixing with the output of a local oscillator 824. The average power of the downcoverted signal is 
monitored by an initial power control 825 that adjusts the average power to the specifications required by the 

45 sampler 826. The adjusted downconverted signal is then sampled in sampler 826 to produce bits of information. 
The bits are then processed In the deprocessor 828, described in more detail in Figure 14b. An estimate of 
the phase differential is taken in the channel estimator 830, as described in more detail in relation to Figure 7 
above, and the estimated phase differential is supplied to a decoder/demodulator 832 to correct the received 
bits. The estimated phase differential Is also supplied to a pre-dlstorter 834 in the transmitter. At the transmitter 

50 in the Base Station, the same blocks are incorporated as in the portable transmitter except that a pre-distorter 
is used to alter the phase of the D8PSK symbols to make the channel appear Gaussian (Ideal) as opposed to 
a fading channel. The initial power control 825 also sends a signal to the pre-distorter 834 to adjust the trans- 
mitted power to an appropriate signal level for the sampler 826 in the first transceiver. It will be appreciated 
that a pre-distorter will be included in the first transceiver's transmitter but that it will not be operable, except 

55 when the first transceiver is operating as a base station. 

Figure 13c shows the functional blocks of the receiver of the first transceiver, which Is the same as the 
receiver in the second transceiver except it does not include an estimator. The processor is illustrated in Figure 
14a and 14c and the deprocessor in Figure 14b and 14c. The processor first Inverse Fourier transforms the 
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128 D8PSK symbols output from the modulator. The transformed symbols are then triplicated as a group so 
that the total number of samples is tripled (see the left side of Figure 4c), with three consecutive groups each 
consisting of the 128 transformed symbols. Next, the three groups are windowed by a Raised Cosine window 
with a roll'Off of 0.25 centered in the middle of the three groups. In other words, the processor takes D8PSK 

5 symbols in, pulse shapes them and inverse Fourier transforms them. On the other hand, the deprocessor un- 
does what the processor did, i.e. it removes the pulse shaping, then Fourier transforms the received signal to 
obtain the original D8PSK symbols. The first two blocks in Figure 14b are similar to the second two blocks in 
Figure 14a except for two differences as follows. In the first block shown in Figure 14b, the repeated groups 
of symbols are partially overlapped, as shown in Figure 14c. In the second block, a rectangular window is used 

10 instead of the Raised Cosine to produce 128 output samples corresponding to the 416 input samples. 

The phase estimator is the same as that shown in Figure 7, except that there are only 1 28 input samples, 
and the same description applies. 

For both the LAN and cellular networks, the present system is designed to operate as a spread spectrum 
system preferably over such bands as are permitted, which at present are the 902 - 928 MHz band, 2.4- 2.4835 

15 GHz and 5.725 - 5.85 MHz. The carrier frequency in the local oscillator shown in Figures 5a, b and c will then 
be 915 MHz in the case of the 902 - 928 MHz band, and the frequencies used for modulation will be centered 
on this carrier frequency. 

Alternative Embodiments 

20 

A person skilled in the art could make immaterial modifications to the inventton described and claimed in 
this patent without departing from the essence of the invention. 

For example, a system may consist of one or more central controllers (comparable to the Base Stations 
in the exemplary celluiar system described) and some slave units (comparable to the portables). The slave 
25 unit executes the commands it receives from the central controller. The cornmands may be requesting the slave 
unit to transmit a receive acknowledge, a status code or information that the slave has access to. The command 
may also be to relay the comnnand or the information to another slave unit 
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#of 
layers 


urban 
area « 
in Km'^ 


indoor 
area. 
in Km^ 


CIR 
in dB 


7 


2 


1.43 


0.02 


13.8 






(C) 










Table IV 





16 




5 



10 



. 15 



20 



25 



30 



35 



40 



45 



50 



CD 



EP 0 562 868 A2 



6.18 Kbps 


9.15 Kbps 


18.77 Kbps 


Vocoder 
Rate 




CO 
G> 


CO 
CJl 


00 


3J O II 

3 ^ 

CD 0,0 


o 
o 


2856 


1995 


1008 


# of FDvc 
BS Tx/Rx 
(21X1XN) 


MHz 


00 




ro 
-f^ 


N=#of 

FDvc/ 

frame 


CJl 

o 

^ . 


2856 


1974 


1008 


#of FDvc 
BS Tx/Rx 
(21X2XN) 


IHz 


00 


CJI 


ro 


N=#of 

FDvc/ 

frame 


5.86 


1008 


CJ> 
CO 

o 


ro 
cn 
ro 


#of FDvc 
BS Tx/Rx 
(7X1 8XN) 


MHz 



17 



EP 0 562 868 A2 




EP 0 562 868 A2 



APPENDIX A 

Copyright § 1991 University Technologies Inc. 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCC.CCCCCCCCCCCCCCCCCCCC 

ccc 

? c 
c c 

C A LISTING FOR THE "PROCESSOR- C 
0 WRITTEN IN FORTRAN C 

C 0 
GCCCCCCCCCCCCCCCCCCCCCCGCCCCGCCCCCCCCCGCCCCCCCCCC 

ccc 

C^^^^--M-^--H^^-+++++++++^^.^-M-+++^H■-l . ^ ^ ^ 
G NAME : BOXNRM(ISTRM) - FUNCTION 

C + 
C FUNCTION : SIMULATES NORMALY DISTRIBUTED DEVIATES 

C ( POLAR METHOD OF BOX & MULLER IS USED) + 

C REF: + 

C ^. • • . 

G GALLING PFiOG / SUBPROG : + 
C + 
G GALLED SUBPROG : URAND + 
C + 

G INPUTS : ISTRM - PASSED FROM THE CALLING ROUTINE + 

C + 

C OUTPUTS : STD. NORMALY DISTRIBUTED RANDOM DEVIATES 

+ ■ .. ■ 

C + 
G PRINCIPAL VARIABLES : + 



G 



+ 



C 

G DATE : 26 APRIL , 1985: + 

G 



+ 



FUNCTION BOXNRM(ISTRM) 

G 

C DECLARATION: 



REAL URAND. VI. V2. R, FAG, GSET. BOXNRM 
INTEGER ISTRM, ISET 
SAVE GSET 



DATA ISET/0/ 

G 

G SHOULD WE GENERATE A PAIR OF DEVIATES '> 
C 

IF(ISET .EQ. 0) THEN 

G 

C YES. WE DONT HAVE AN EXTRA DEVIATE HANDY. SO PICK 
G TWO UNIFORM NUMBERS IN THE SQUARE EXTENDING FROM 
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C -1 TO -t-1 IN EACH DIRECTION 

C 

1 VI = 2. • URAND(ISTRM) - 1 . 

V2 = 2. * URAND(ISTRM) - 1 . 

0 

C SEE IF THEY ARE IN THE UNIT CIRCLE, IF THEY ARE NOT GO TO 1 
C 

R = Vf2 + V2"2 
IF(R.GE. 1.)G0T01 

C 

0 NOW MAKE THE BOX-MULLER TRANSFORMATION TO GET TWO 

NORMAL DEVIATES 

C 

FAG = SQRT(-2. * ALOG(R)/R) 

C 

C ' RETURN ONE AND SAVE ONE FOR NEXT TIME 
C 

GSET = VI • FAC 
BOXNRM = V2 * FAC 

C 

C SET FLAG 
C 

ISET = 1 

C 

ELSE 

C 

C WE HAVE AN EXTRA DEVIATE HANDY. SO RETURN IT 
C 

BOXNRM = GSET 

C 

C UNSET THE FLAG 
C 

ISET = 0 

C 

ENDIF 

0 

RETURN 
END 
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subroutine clemod16(cielta.cc) 
complex cc 

common /pie1/piover4.piover8,piover1 6 

lf(delta.ge.31.*piover16.or .delta.lt. 1.'piover16)cc=(+1.+1) 
if(delta.ge. 1 /ploverl 6.and.delta.lt. 3.*piover1 6)cc=(+1 ,-l ) 
if(delta.gG. 3.*piover1 6.and.delta.lt. 5.*plover16)cc=.(-l,-l) 
if(delta.ge. 5/piover1 6.and.delta.lt. 7.*pioverl6)cc-(-1,+l) 

if{delta.ge. 7.*piover16.and.delta.lt. 9.*pioverl6)cc=(-1 ,+3) 
if(delta.ge. 9.*plover16.and.delta.lt.1 1.*piover16)cc=(-1.-3) 
if(deita.ge.1 1 .•piover16.and.delta.lt.13.''piover16)cc=(+1 ,-3) 
if(delta.ge.13.*piover16.and.deIta.lt.15.*ploven6)cc=:(+1 ,+3) 

if(delta.ge.15.*piover16.and.delta.lt.17.*piover16)cc=(+3.+3) 
if(delta.ge.17.*piover16.and.delta.lt.19.*piov6r16)cc=(+3,-3) 
if(delta.ge.19.*piover16.and.delta.lt.21.*piover16)cc=(-3.-3) 
if(delta.ge.21 .*piover1 6.and.delta.lt.23.*plover1 6)cc=(-3.+3) 

if(delta.ge.23.*piover1 6.and.delta.lt.25.*piover1 6)cc=(-3.+1 ) 
if(delta.ge.25.'piover1 6.and.delta.lt.27.*piover1 6)cc=(+3,+1 ) 
if(delta.ge.27.*piover1 6.and.deIta.lt.29.*piover1 6)cc-(+3,-1 ) 
if(delta.ge.29.*piover1 6.and.delta.lt.31 .'pioverl 6)cc=(-3,-1 ) 

return 
end 
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subroutine demod4(delta,cc). 
complex cc 

common /pie1/piover4,piover8,piover1 6 

if(delta.ge.0.*piover4.and.delta.lt.2.*piover4)cc=cmplx(+1 +1) 
if(delta.ge.2.*piover4.and.delta.lt.4.*piover4)cc=cmplx(+l'-i) 
if(delta.gG.4.*piover4.and.delta.lt.6.*piover4)cc=cmplx(-l,-i) 
if(delta.ge.6.*piover4.and.delta.lt.8.*piover4)cc=cmplx(-1,+l) 

return 
end 
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subroutine demod8(delta,cc) 
complex CO 

10 common /pie1/piover4.piover8.piover16 

if(deltageJ5/piov^ 1.*P»over8)cc=(-h1 +1) 

if(de ta.ge. l/pioverS.and.deltaJt. 3/piover8)cc=(+lVi 
•IS^ I^-^®- 3^;piover8.and.deltaJt. 5/piover8)cc=H .-i ) 
^5 if(delta.ge. 5/piover8. and.delta.lt. 7.*piover8)cc={.l i+i ) 

ilSfllf'?!* 9.*piover8)cc=(-3,+1) 

20 




return 
end 
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C SUBROUTINE FFT{) - REFER TO THE BOOK BY 'GONZALEZ & WINTZ* 
c 

c f is the input array to the fft 

c f is also the output array from the fft with +ve freq. at the begining 
^ and -ve freq. at the end. 

c f1 is the output array from the fft with -ve freq. at the begining 
^ and +ve freq. at the end. 

c ind is the index array for f 1 . 

0 if ix=0 the F is the input array, 
c if ix=1 the F1 is the input array. 

subroutine fft (ix.f.lnd.fl , In) 

integer In, n, nv2, nm1, j, i. k. lie. Ie1. ip.ind(2048) 

real pi ' 

complex f(2048). f1(2048).u. w, t 

c 

pl = 4.0 *atan(1.0) 
n = 2"ln 
nv2=.n/2 
nm1 = n-1 

if(lx,eq.1)call deswap(f,f1,n) 
j = 1 

do3i=1.nm1 

(i ge j) go to 1 
t = f(j) 
m = f(i) 
f(i) = t 

1 k = nv2 

2 if (k.ge.j) go to 3 

j=i-k 

k = k/2 
go to 2 

3 j = J+k 

do 51=1, In 
le = 2**1 
lei = le/2 
u = (1.0, 0.0) 

w = cmplx( cos(pi/!el ), -sin(pi/le1 ) V 
do5j=1,le1 
do4i=j, n, le 
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ip = i + Ie1 
t = f(ip)*u 
f(ip) = f(i)-t 

4 f(i) = f(i)+t 

5 u = u * w 
do 6 i= 1. n 

f(i) = f(i) / float (n) 
15 6 ind(i) = i-n/2-1 

call swap(f.f1,n) 

return 

end 
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5 



10 



c 
c 
c 



Inverse fft transform 



c 



c if ix«0. f is the input array to the Ifft with +ve freq, at the begining 



c if ix=1 , f 1 is the input array to the Ifft with -ve freq. at the begining 
c and +ve freq, at the end. 

subroutine tdft (ix,xft.ind,xft1 .nft) 
integer nft, i, mft,lnd(2048) 

complex xft(nft).xft1 (2048), buff(2048), buffi (2048) 



mft = nint(alog10(fioat(nft))/alog 10(2.0)) 
print *, • nft and mft are:*, nft, mft 

if(ix.eq.1 )call deswap(xft,xft1 .nft) 

, do 7 i= 1 . nft 

7 buff(i) = conjg(xft(i)) 

c ffl() does 1-D Fourier transform 

call fft(0. buff. ind. buffi, mft) 

do 8 i= 1 , nft 
xft(l) = nft*conjg(buff(i)) 

8 xft1 (i) = nft*conjg(buff 1 (i)) 



15 



c 



and -ve freq. at the end. 



c 



40 



return 
end 
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subroutine swap(f,f1,n) 
complex f(2048). f 1(2048) 

do i=-n/2r1 

f1(i+n/2+1) =f(n+i+1) 
enddo 

do i=0,n/2-1 

f1{i+n/2+1) = f(i+1) 
enddo 

return 
end 



subroutine deswap(f»f1,n) 
complex f(2048). f 1(2048) 

do i==-n/2,-1 

f(n+i+1) = f1(i+n/2+1) 
enddo 

do i=0,n/2-1 

f(i+1) = f1(i+n/2+1) 
enddo 

return 
end 
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real pi,TWOPI,arRe(0:51 2).arlm(0:51 2),norm.p. 

10 ' * curphase(0:512).oldph.ph.psi(0:512).RC(-2048:2048) 
integer ind(2048),counter 

complex 2(0:51 2).r(0:512).c(0:512).d(0:512),cc 
complex F(2048),F1 (2048) 
character strlng'50 
IS common /pie/pi 

common /pie1/piover4,plover8,pjoverl6 
common /dec/no 

open(1 00,file='transmitted') 
20 open(200,file='received') 
open(400.file='BER') 
open(800.file=*envelope') 
open(310,file=:'RC') 
open(350.file=*data') 
open(351 ,file='random') 
pi = 4. * atan(1.G) 
piover4 = pi/4. 
piover8 = pi/8. 
pioverlS = pi/16. 
TWOPI = 2. * pi 
read(*,102)string 
write(*,102)strlng 
read(*,'')number 
35 readC.*)M 
readC/)MM 
read(*,')no 
102 format(a50) 
C open(101,file=string) 



25 



30 



40 



numb=255 
LL= MM/2 -1 
counter = 0 

call raisedcos(MM,RC) 



45 

C do 1=1 .numb 



c***"" READING 

C read(101,*)p.arRe(i).arlm(i) 

C norm = ( 900. / p )**2 

C arRe(i) = arRe(i)/norm 

C arlm(i) = arlm(l)/norm 

C z(l)=cmplx(arRe(j),arlm(i)) 

C write(800.*)i.sqrt(real(Z(i)*conjg(Z(i)))) 

55 C enddo 

C z(0) =cmplx(arRe(1),arlm(1)) 

e z(MM)=cmplx(arRe(numb).arlm(numb)) 
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do ii=1, number 
C""" FIRST LOOP 

curphase(O) = urand(3)*TWOPI 

d(0) = cmplx(cos(curphase(0)).sin(curphase(0))) 

C* STARTING THE DATA TRANSMISSION 
C"*"» DATA GENERATION 
do 10 1=1, MM, 1 

if(M.eq.4)call modu4(cc.delta) 

if(M.eq.8)call modu8(cc,delta) 

if(M.eq.1 6)call modul 6(cc,delta) 

write(500,*)l.cc.delta 

c(l) = cc 

psi(l) = urand(4)'TW0PI 
write(351,*)psi(l) 

curphase(l) = curphase(l-1 ) + delta + psi(l) 
d(l) = cmplx(cos(curphase(l)).sin(curphase{l))) 
enddo 

C""" MODULATOR 

do k=0.MM,1 

write(1 00.*)k,sqrt(real(d(k)*conig(d(k)))) 
0 r(k)-d(k)*2(k) 
r(k) = d(k) 
write(451 .•)k,r(k) 
enddo 

C FIND THE TIME RESPONSE * 
j=0 

dok=1.MM,l 
j=j+1 

F10)=r(k-1) 
enddo 

CALL tdft(1 ,F.lnd,F^1 ,MM) 

open(300,file='time_RC5SP) 

open(301,flle=»'time0') 

open(302.file='time_NOSP) 

call wlndow(MM,LL,RC,Fl, counter) 

C gotolOI 

C DEMODULATOR 

oldph = atan2(aimag(r(0)),real(r(0))) 
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dok=1,MM-1 

ph = atan2{aimag(r(k));real(r(k))) 

delta = ph - oldph - psi(k) 

if(delta.gt.TWOPI)delta=delta-TWOPI 

if(delta.gt.TWOPI)delta=delta-TWOPI 

lf(delta.gt.TWOPI)delta=delta-TWOPI 

if(delta.lt.0.)delta=delta+TWOPI 

if(delta.lt.0.)delta=delta+TWOPI 

if(delta.lt.0.)delta=delta+TWOPI 

if(M.eq.4)call demod4(delta.cc) 
if(M.eq.8)call demod8(delta,cc) 
if(M.eq.l6)call demodl6(delta.cc) 
write{501,*)k.cc,delta 

if(real(cc).ne.real(c(k)))errorRe=errorRe+1 
lf(aimag(cc).ne.aimag(c(k)))errorlm=errorlm+1 

write(400,')k,* •,real(cc),real(c(k)), 
& • ',aimag(cc),aimag(c(k)) 

oldph = ph 
enddo 

prinf.'errorRe=*.errorRe.' '.•errorlm^'.errorlm 
10 continue 
101 continue 

enddo 

rewind(301) 

call quantize(LL,counter.prod) 
write(306,*)prod 



stop 
end 
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10 



15 



20 



25 



subroutine modu16(cc,delta) 
complex cc 

common /pie1/piover4,piover8,piover1 6 

si = urand(1) 

s2 = urand(2) 

S3 = urand(l) 

s4 = urand(2) 

if(s1.ge.0,5)s1==-h1 

if(s1.lt.0.5)s1=-1 

If(s2.ge.0.5)s2=+1 

if(s2.lt.0.5)s2=-1 

if{s3.ge.0.5)s3=+1 

If(s3.lt.0.5)s3=-1 

if(s4.ge.0.5)s4=-h1 

if(s4.lt.0.5)s4=-1 

cc = cmplx(s1+2*s3,s2+2*s4) 

write(350.*)cc * 



30 



35 



40 



45 



if(cc.eq. 
if(cc.eq. 
if(cc.eq. 
if(cc,eq. 

if(cc.eq. 
if(cc.eq. 
if(cc.eq. 
if(cc.eq. 

if(cc.eq. 
if(cc.eq, 
if(cc.eq. 
if(cc.eq, 

if(cc.eq. 
if(cc.eq. 
if{cc.eq. 
if(cc.eq. 



+1 ,+1 ))delta=0/piover8 
+1 ,-1 ))delta=1 .*piover8 
-1 ,-1 ))delta=:2/piover8 
-1,+1))delta=3/piover8 

1 .+3))delta=4.*piover8 
-1 .-3))delta=5/piover8 
+1 r3))delta=6.*piover8 
+1 ,+3))delta=7.*piover8 

-h3,+3))delta=8.*piover8 
4-3.-3))delta=9.*piover8 
-3.-3))delta=10.*piover8 
-3,+3))delta=1 1 .*piover8 

-3,+1 ))delta=1 2.*piover8 
+3,+1 ))delta=1 3.*piover8 
+3,-1 ))delta=1 4.*piover8 
3,-1))delta=15/piover8 



50 



return 
end 
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subroutine modu4(cc,delta) 
complex cc 

common /pie1/piover4,piover8,piover1 6 

s1 = urancl(l) 
s2 = urand(2) 
if(s1,ge.0.5)s1=+1 
if(s1.lt.0.5)s1=-1 
if(s2.ge.0.5)s2=+1 
if(s2.lt.0.5)s2=-1 
cc = cmplx(s1,s2) 



lf(cc.eq.(+1 .+1 ))delta=1 /piover4 
if(cc,eq.(+1 .-1 ))delta=3/piover4 
if{cc,eq,(-1 .-1 ))delta«5.*piover4 
if(cc-eq.(-1 .+1 ))delta=7/piover4 

return 
end 
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write(350/)cc 
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subroutine modu8(cc,delta) 
complex cc 

common /pie1/piover4,piover8,piover1 6 

s1 = urand{1) 

s2 = urand(2) 

s3 = urand(1) 

if(s1.ge.0.5)s1=+1 

if(s1.lt.0.5)s1=-1 

if(s2.ge.0.5)s2=+1 

If(s2.lt.0.5)s2=-1 

if(s3.ge.0.5)s3=+1 

if(s3.lt.0.5)s3=-1 

cc = cmplx{s1+2*s3,s2) 

pp = urand(4)*TWOPI 

write(350,*)cc 

if(cc.eq.(+1 ,+1 ))delta=0.*piover4 
if(cc.eq.(+1 .-1 ))delta=1 .*piover4 
if(cc.eq.(-1 .-1 ))delta=2.*piover4 
if(cc.eq.(-1 .+1 ))delta=3/piover4 

if(cc.eq.(-3.+1))delta=4.*piover4 
if(cc.eq.(+3,+1))delta=5/piover4 
if(cc.eq.(+3,-1))delta=6.*piover4 
if(cc.eq.(-3,-1))delta=7.*piover4 

return 
end 
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subroutine quantize(LL,counter,prod) 
integer Fre,Fim,counter 



xmax=-100. 
xmin=100. 
do j=1 , counter 
15 read(301.*)i.arRe,arlm 

if(arRe.gt.xmax)xmax=arRe 
if(arlm,gtxmax)xmax=arlni 
if(arRe.ltxmin)xmin=arRe 
if(arlmJt-xnnin)xmin=arlm 
20 enddo 

if(abs(xmax).gt.abs(xmin))xxx=abs(xmax) 
if(abs(xmax).lt.abs(xmln))xxx=abs(xmin) 
prod = 1 28./XXX 

sumRe = 0. 
sumlm = 0. 
surRe = 0. 
surim = 0. 
rewind(301) 
do j=1. counter 
read(301 /)i,arRe,arlm 
Fre = 1 28 + prod*arRe 
Hm = 128 + prod*arlm 
35 C if(Fre.gt.255)Fre-255 

C if(Fre.lt.0)Fre=0 
C if{Rm.gt.255)Fim=255 
C if(nm.lt.O)Fim=0 

write(40l/)i,Fre,Fim 
40 C wrjte(402.*)j.Fre 

C write(403,*)j,Fim 

sumRe = sumRe + (arRe-(Fre-1 28.)/prod)**2 
sumlm = sumlm + (arlm-(Flm-128.)/prod)**2 
surRe = surRe + (arRe)**2 
surIm = surIm + (arlm)**2 
enddo 

print*/quanti2ation noise=, (sumRe + sumlm)/{surRe + surIm) 
print*,prod ' 



return 
end 
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subroutine raisedcos(MM,RC) 
real RC(-2048:2048).pi 
common /pie/pi 

do i=-MM/2-MM/8.-MM/2+MM/8 
„ t = -float(i)/float(MM) 

RC(i)=(0.5*(1.-sin(4/pi*(t-0.5)))) 
write(310.*)i.RC(i) 
enddo 

20 do i=-MM/2+MM/8+1.MM/2-MM/8-1 

t = float(i)/float(MM) 
RC(i)=1. 

write(310.*)i.RC(i) 
enddo 

25 

do NMM/2-MM/8.MM/2+MM/8 
t = float(i)/float(MM) 
RC(i)=(0.5*(1.-sin(4.*pi*(t-0.5)))) 
write(310.*)i.RC(i) 
^ enddo 

return 
end 

35 
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10 subroutine demod4(delta.wRe.wlm) 

y = cmplx{cos(delta),sin(delta)) 
lf(real{y).ge.O)then 

wRe=+1. 
IS else 

wRe=-1. 
endif 

if(aimag{y),ge.O)then 
2^ wlm«-h1. 
else 

wlm=-1. 
endif 

return 
end 
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C NAME : URAND(ISTRM) - FUNCTION 

^ + 

C FUNCTION : GENERATES STANDARD UNIFORMLY DISTRIBUTED 

C RANDOM NUMBERS . + 

C ( USES THE RECURSION: SEED(ISTRM) = + 

^ 1 6807 • SEED(ISTRM) MOD (2"(31 ) - 1 ) ) + 

J SOME COMPILERS REQUIRE THE DECLARATION: 

C INTEGERS ISTRM. K1 + 

C REF: V 

C CALUNGPROG/SUBPROG : ^ 

C + 

C CALLED SUBPROG :NIL . 

C + 

C INPUTS : ISTRM - THE STREAM NUMBER + 

^ 0<SEED(ISTRM)< 2147483647 + 

^ + 

C OUTPUTS : STD. UNIFORMLY DISTRIBUTED RANDOM NUMBERS 

^ NEW VALUE OF SEED(ISTRM) + 

C PRINCIPAL VARIABLES : 

C DATE : 26 APRIL. 1985. + 

c 




^ FUNCTION URAND{ISTRM) 

C DECLARATION : 
C 

INTEGER SIZE3 
PARAMETER(SIZE3 = 20) 
REAL URAND 

INTEGER SEED(SIZE3).ISTRM.ISTRG.K1 
save seed 

data seed(1)/456789/ 
data seed(2)/l 037625857/ 
data seed(3)/1 02899405/ 
data seed(4)/535417/ 
data seed(5)/21 1 229747/ 
data seed(6)/2952853/ 
data seed(7)/2987570/ 
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data seed(8)/45376256/ 
data seed(9)/1 7890857/ 

C 

" C GENERATE A U(0.1 ) VARIATE 



ISTRG = ISTRM 

K1 =SEED(ISTRG)/1 27773 

SEED(ISTRG) - 16807 * ( SEED(ISTRG) - K1 • 127773) - K1 * 2836 
IF( SEED(ISTRG) LT. 0) SEED(ISTRG) = SEED(ISTRG) + 2147483647 
URAND = FLpAT(SEED(ISTRG)) • 4.65661 2875E-10 

RETURN 
END 
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subroutine window(MM.LL.RC.F.counter) 

complex F(2048).FO(-2048:2048).Fm(-2048:2048).Fp(-2048-2048) 
Fsum(-2048:2048).FRC(-2048:2048) '* 
real RC(-2048:2048) 
integer counter.counti 
common /dec/no 
save counti 
data countl/0/ 

counti = counti + 1 " 
do i=-M!VI-LL,MM+LL+1 

FO(i) = (0..0.) 

Fm(i) = (0..0.) 

Fp(i) = (0..0.) 
enddo 

do i=-LL,LL+1 

FO(i) = F(i+LL+1) 

Fm(i-MM) = FO(i) 

Fp(i+MM) = FO(i) 

write(304/)i.F0(i) 
enddo 

do i=-MM-LL,MM+LL+1 

Fsum(i) = FO{i) + Fm(i) + Fp{i) 

^^wnte(302.*)i.sqrt{(real(Fsum(i)))-2+(aimag(Fsum(i)))"2) 

do i=-LL - {LL+1)/4-8.-LL - (LUl)/4-1 
counter = counter + i 
FRC(i) = cmplx(0.,0.) 

write(301 ,*)counter,0.,0. 
enddo 

do i=-LL - (LL+1 )/4.LL+1 + (LL+1 )/4 
counter = counter + i 
FRC(i) = Fsum(l)*RC(i) 

write(300*)counter.sqrt((real(FRC(i)))"2+(aimag{FRC(i»)**2) 

wnte(301.*)counter.real(FRC(i)).aimag(FRC(i)) 

write(202.*)counter.real(FRC(i)) 

write{203,*)counter,aimag(FRC(i)) 
enddo 

do i=LL+1 + (LL+1 )/4+1. LL+1 + {LL+1)/4+8 

counter = counter + i 

FRC(i) = cmplx(0..0.) 

write(301 .*)counter.O..O. 
enddo 
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if(count1.eq.no)then 
,0 do i=-LL - (LL+1)/4.8.LL+1 + (LUl)/4+8 

counter = counter +1 
FRC(i) = cmplx(0..0.) 
write(301 ,*)counter,0.,0. 
write(202,*)counter,real(FRC(i)) 
15 write(203,')counter.aimag(FRC(i)) 
enddo 
endif 



return 
end 
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FFU\GS = 

FFILES = main.! urand.f fft.f window.! raisedcos.f quantize.! \ 
niodu4.f moudS.f modulG.f demod4.f demodS.f demodlG.f 

OFILES = main.o urand.o fft.o window.o raisedcos.o quantize.o \ 
modu4.o moduB.o modu16.o demod4.o demodS.o demod16.o 

a.out: ${OFILES} 

f77 ${FFLAGS} ${OFILES} 
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.7..Ay../../HF_LF/TWOCHANNELy900/test.1 

9 

16 

64 

9 
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mTtntrtru.mu»irin tfrnmmrttfri it.t.r,r,rr- 

C WRITTEN IN FORTRAN C ^ 

C c 
C A LISTING FOR THE "DEPROCESSOR" C 
^ C 

S$5£^^^^^^^^^^^^cccccccccccccccccccccccccccccccc 

subroutine demod16(delta,cc) 
complex cc 

common /piel/piover4,pi6ver8,pioverl6 

lf{delta.ge.31 .*piover1 6.or .delta.lt. 1 .*plover1 6)cc={+1 +1 ) 
if(delta.ge. l.*piover16.and.delta.lt. 3.*piover16)cc=(+1,-1) 
lf(delta.ge. 3.*piover16.and.delta.lt. 5.*plover16)cc=(-1 -1) 
if(delta.ge. 5.*piover1 S.and.delta.lt. 7.*piover1 6)cc={-1 [+1 ) 

if(delta.ge. 7.*piover1 6.and.delta.lt. 9.*piover16)cc=(-1 +3) 
if(delta.ge. 9.*piover1 s.and.delta.lt. 1 1 .'pioverl 6)cc=(-1 -3) 
if(delta.ge.1 1 .'pioverl 6.and.della.lt.1 3.*piover1 6)cc=(+l -3) 
if(delta.ge.1 3.*piover1 S.and.delta.lt.l 5.*piover1 6)cc=(+1 ,'+3) 

if(delta.ge.1 5.*ploverl e.and.delta.lt.l 7.*piover1 6)cc=(+3 +3) 
if{deita.ge.1 7.*piover1 e.and.delta.lt.l 9.'piover1 6)cc=(+3'-3) 
if(delta.ge.1 9.*piover1 6.and.delta.lt:21 .*plover1 6)cc=(-3 -3) 
if(delta.ge.21.*ploverl6.and.delta.lt.23.*pioverl6)cc=(-3.+3) 

if(delta.ge.23.*piover16.and.delta.lt.25.*piover16)cc=(-3 +1) 
if(delta.ge.25.*piover1 6.and.delta.lt.27.*piover1 6)cc=(+3 +1 ) 
if(delta.ge.27.*piover16.and.delta.lt.29.*pioverl6)cc=(+3'-1) 
if(delta.ge.29.*pioverl6.and.delta.lt.31.*piover16)cc=(-3.-1) 

return 
end 
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subroutine demod4((lelta,cc) 
complex cc 

common /piel /piover4,piover8,piover1 6 

if{delta.ge.0.*piover4.and.delta.lt.2/piover4)cc=cmplx(+1 +1) 
if{delta.ge.2/piover4.and.delta.lt.4/piover4)cc=cmplx(+l '-1 ) 
if(delta.ge.4.*piover4.and.delta.lt.6.*piover4)cc-cmplx(-1 -1 ) 
if{delta.ge,6/piover4.and-delta.lt.8/pjover4)cc=cmplx(-l i+1 ) 

return 
end 
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subroutine demod8(delta,cc) 
complex cc 

common /pie1/piover4,piover8,piover1 6 

if(delta.geJ 5/piover8.or .delta.lt. 1 /piover8)cc=(+1 .+1) 
if(delta.ge. 1 /pioverS.and.delta.lt. 3.*piover8)cc=(+1 ,-1) 
if(delta.ge. 3.*piover8.and.delta.lt. 5.*piover8)cc={-1,-1) 
if(delta.ge. 5.*piover8.arid.delta.lt. 7."piover8)cc=(-1,+l) 

if(delta.ge. 7.*piover8.and.deltaJt. 9.*piover8)cc=(-3.+1) 
if(della.ge. 9.*piover8.and.delta.lt.1 1.*piover8)cc=(^-3,+1) 
if(delta.ge.1 1 .•plover8.and.delta.lt. 1 3.*piover8)cc=(+3.-1 ) 
if(delta.geJ3.*piover8.and.deltaJt.iS.*piover8)cc=(-3.-1) 

return 
end 
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subroutine dewindow{MM,numb,F) 
integer shift 

complex F(2048).FO(-2048:2048).Fm(-2048:2048). 
& Fsum{-2048:2048).Fp(-2048:2048).fprime 
^« common /SHHH/shift 

common /SDHH/prod 

do 1=1, numb 
read(l 01 ,*)j,arRe.arlm 
20 F(i) « cmplx(arRe-shift.arlm-shift) 

enddo 

if(coun1.eq.no)then 
do i=1.numb 
25 read(101,*)j,arRe,arlm 
enddo 
endif 

LUMM/2-1 

30 do j=-MM-numb/2+1 ,MM+numb/2 

FO(i) = (0.,0.) 
Fm(i) = (O.A) 
Fp(i) = (0.,0.) 
enddo 

35 

do i=-numb/2+1, numb/2 
FO(i) = F{i+numb/2) 
Fm(i-MM) = FO(i) 
Fp(i+MM) = FO(i) 
40 enddo 

do i=-MM-numb/2+1.MM+numb/2 

Fsum{i) « FO(i) + Fm{i) + Fp{i) 
enddo 

45 

sum = 0, 

sup = 0. 

do i=-LL.LL+1 

write(303/)i.sqi1({real(Fsum(i)))**2 
50 & +(aimag{Fsum(i)))**2) 

F(i+LL+1) = Fsum(i) 

write(304/)i.Fsum(i) 

C read(604/)j,Fprime 
55 C sum = sum + conjg(Fsum(i)/prod-Fprime)* 

C & (Fsum(i)/prod-Fprime) 
C write(602,*)i,abs(real{Fsum(i)/prod)) 
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C & -abs(real(Fprime)) 

C write(603/)i,abs(aimag(Fsum(i)/prod)) 

C & -abs(aimag(Fprime)) 

C write(605,*)i,Fsum(i)/prod-Fprime 

C 

C sup = sup 4- conjg(Fprime)*(Fprime) 
enddo * 
20 C print* /quantization error= ,sum,sup,sum/sup 

return 
end 
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c SUBROUTINE FFT() - REFER TO THE BOOK BY 'GONZALEZ & WINTZ 
c 

c f is the input array to the fft 

c f is also the output array from the fft with +ve freq. at the begining 
c and -ve freq. at the end. 

c f1 is the output array from the fft with -ve freq. at the begining 
c and +ve freq. at the end. 

c ind is the index array for f 1 . 

c if ix=0 the F is the input array, 
c if ix=1 the F1 is the input array. 



subroutine fft (ix.f.ind.f 1 . In) 

integer In. n, nv2, nml, j. i. k. le. let. ip.ind(2048) 

real pi 

complex f(2048), f 1 {2048),u, w. t 

c 

pi - 4.0 * atan{1 .0) 
n = 2**ln 

nv2 = n/2 
nml = n-1 

if(ix.eq.1)call deswap(f.fl.n) 
j=1 

do 3 i= 1, nml 
if (i.ge.j) go to 1 
t = fa) 
fG) = f(i) 
f(i)=t 

1 k = nv2 

2 if (k.ge.j) go to 3 
j = j-k 

k = k/2 
goto 2 

3 j = j+k 

do 5 1= 1, In 
le = 2**1 
lei = le/2 
u = (1.0. 0.0) 

w = cmplx( cos(pi/le1), -sin(pi/lGl) ) 
do5j=l.le1 
do 4 i= i. n. le 
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ip = i + lei 
t = f(ip) • u 
f(ip) = f(i)-t 
f{i)=f(i) + t 

u = u • w 

do 6 i= 1 , n 

f(i) = f(i) / float (n) 

ind(i) = l-n/2-1 

call swap(f,fl,n) 

return 

end 
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C" • * 

C Inverse fft transform 

c 

c if ix=0, f is the input array to the Ifft with +ve freq. at the begining 
c and -ve freq. at the end. 

c if ix=1 , f 1 is the input array to the Ifft with -ve freq. at the begining 
c and +ve freq. at the end. 

subroutine tdft (ix.xft.ind.xftl ,nft) 
integer nft, I, mft.ind(2048) 

complex xft(nft).xft1 (2048). buff(2048), buffi (2048) 

c 

mft = nint(alog10(float(nft))/alog1 0(2.0)) 
print *, ' nft and mft are:', nft, mft 

if(ix.eq.1 )call deswap(xft,xft1 .nft) 

do7i=1,nft 

7 buff(i) - CGnjg(xft(i)) 

c fft() does 1 -D Fourier transform 

call fft(0.buff,ind.buff 1 , mft) 

do 8 i= 1 , nft 
xft(i) = nft*conjg(buff(i)) 

8 xft1 (i) = nft'conjg(buff1 (i)) 

return 
end 
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subroutine swap(f,fl ,n) 
complex f(2048). f1 (2048) 

do l=-n/2,-1 

f1(i+n/2-».1)-f(n-»-i+1) 
enddo 

do i=0.n/2-1 

f1(i+n/2+1) = f(l+1) 
enddo 

return 
end 
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subroutine deswap(f,fl,n) 
complex f(2048). f1 (2048) 

do i=-n/2,-1 

f(n-i-i+1) = f1(i+n/2+1) 
enddo 

do i=0.n/2-1 

f(i+1) =f1(i+n/2+1) 
enddo 

return 
end 
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real pi.TWOPI,oldph,ph,psi(0:512) 

complex r(0:512).c{0:512),cc,cl 

complex F(2048),F1(2048) 

integer ind(2048), shift 

character stringMO.stringI M0.string2*40 

common /pie/pi 

common /piel /piover4,piover8,piover1 6 
common /SHHH/shift 
common /SDHH/prod 

open{1 00,file= transmitted*) 
open(200,file=«Yeceived') 
open(400.file=BER') 
open(800.file=*envelope') 
open{310,file=RC') 
pi = 4. •atan(1,0) 
piover4 « pi/4. 
pioverS = pi/8, 
pioverl 6 = pi/1 6. 
TWOPI = 2. ' pi 
read(M02)string 
writer. 102)string 
read(M02)string1 
write{*,102)string1 
read{\102)string2 
write{*,102)string2 
read(*,*)number 
read{*,*)M 
read(*.*)MM 
read{*.*)ni 
102 format(a40) 

open(101,file=string) 
open(350,file=string1) 
open(351 ,file=string2) 
open(604.filev../fort.304') 
open(704.file=»\./fort.501 ') 
shift = 0 
prod = 1 . 

if(string.eq/../fort.40r)shift=128 
if(string.eq/../fort.401')read(306,*)prod 



numb=MM*1.25 + 16 
c READING ' 
sum = 0. 
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do number 

call clewindow(MM,numb,F) 

call m(O.F.ind.Fl.ni) 

i=0 

dok=0.MM-1.1 

r(k) = F1G) 
write(451 .*)k.r(k) 
read(350.')c(k+1) 
read(351,*)psl{k+1) 
enddo 

C—"- DEMODULATOR 

oldph = atan2(aimag(r(0)),real(r(0))) 
dok=1,MM-1 

ph = atan2(aimag(r(k)),real(r(k))) 

delta = ph - oldph - psi(k) + pi 

if(delta.gt.7WOPI)delta=delta-TWOPI 

if(delta.gt.tWOPI)delta=delta-TWOPI 

lf(delta.gt.TWOPI)delta=delta-TWOPI 

if(delta.lt.0.)delta=della+TWOPI 

lf(delta.lt.0.)delta=delta+TWOPI 

if(delta.lt.0.)delta=delta+TWOPI 

if(M.eq.4)call demod4(delta,cc) 
if(M.eq.8)call demod8(delta,cc) 
if(M.eq.1 6)call demodi 6(delta,cc) 
write(501 ,*)k.cc,delta 
read(704,*)l,cl.deltap 
diff » delta-deltap 
if{diff.gt.+pi)diff=diff-TWOPI 
if(diff.lt.-pi)diff=diff+TWOPI 

lf(k.ne.MM/2.and.k.ne.MM/2+1)sum = sum + dlff"2 
write(606.*)k,diff 

lf(real(cc).ne.real(c(k)))errorRe=eiTorRe+1 
if(aimag(cc).ne.aimag(c(k)))errorlm=errorlm+1 

write(400,*)k.* *.reai(cc).real(c{k)). 
& • '.aimag(cc).aimag{c{k)) 

oldph s ph 
enddo , 

prinf.'errorRe=.errorRe.* '.'errorlm^.errorlm 
print*.'sum='.sum/flcat(MM-3) 
10 continue 
enddo 
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Stop 
end 
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FFLAGS - 

. FFILES = main.f fft.f dewindow.f demod4.f demod8.f demodlS.f 

OFILES = main.o fft.o dewindow.o demod4.o demodS.o demod16.o 

a.out: ${OFILES) 

f77 ${FFLAGS} ${OFILES} 
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../fort.401 
../data 
../random 
9 

16 
64 
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SOFTWARE FOR SIMULATION ON DSP32C 

The following software program was developed aiid run using D3SIM of AT&T It 
simulates the operation of Uic equalizer on a DSP32C 



IS 



^include <matliJi> 
^include <stdio.h> 
^Mncludc <libap.h> 



20 
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SO 



4) =( -1 . -l. 



main program 



double twopi«6.283 185372; 
/* List the original transmillcd symbols */ 
intjj2[5 

1 . 

1 , 

1 . 

1 . 

-1 . 

-1 . 

1 . 

-1 . 

1 . • 

-1 . 

1 . 

-1 . 

-1 . 

-I . 

-1 . 

-1 . 

-1 . 

1 . 

-I . 

1 . - 

-1 . 

1 . - 

1 . ■ 

-1 . 

-1 . 

1 . - 

1 , ■ 

-1 . ■ 

1 . - 

-1 . ■ 

-1 . ■ 

-1 . 
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1 . 
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-1 . 
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-1 . 
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1 . 
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-1 . 
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•1 . 
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1 . 


-1. 


1 . 


1. 


1 , 
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1 . 


1. 


-1 . 


1. 


1 . 


1. 


1 . 


1. 
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. 1. 


-1 


. -1. 


-1 


. -1, 
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. -1. 
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. 1. 
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. -1. 
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. -1. 
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. -1. 
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. 1. 
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. -1. 
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. -1. 
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. -1. 
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. 1. 
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. 1. 


-1 


. 1. 


1 


. -1. 


-1 


. 1, 


1 , 


. J. 


-1 , 


. -1. 


I . 


1. 


1 . 


1. 


-1 . 


-K 


-1 . 


1. 




-1, 


-I . 


-1. 


1 . 


1. 


-1 . 


-1, 
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1 , 


1. 


1 , 


-K 


1 . 


.1 

*• 


1 

* > 


.1 
*» 


.1 

* ff 
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•A, 




1 


1 

* t 


if 


1 

* t 


-i» 


1 

A > 


-it 


1 

"A , 


If 


1 
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-1. 


-1 . 


1. 


-1 . 


1. 


-1 , 


1, 


-1 . 


1. 


-1 . 


-1. 


-1 . 
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1 . 


1, 


1 . ■ 
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1 . - 


1. 
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1. 


-1 . 


1. 


-1 . 


1. 


-I . 
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-1. 




-1. 


-1 . 


1. 




I. 


-1 . 
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-1. 
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1. 
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' -1 . -1. 
1 . 1. 
1 . -1. 
1.1. 
1 . -1. 
-1 . -1. 
1 . -1. 
-1 . -1. 

1 . 1 

); 

/* list ihc received I and Q */ 
double 2[514J=( 

1-8826647996902 .0.57904273271561 

-lJ23iSlmfi??° .1.9248703718185. 

nT?I^P"°18967 , 1.1806204319000 

?-i^!?5002686977 . -1.8899723291^7 
-1.3761224746704 .1.2521643638611 
-T 3SllSS?r1??i'^-°2 • -l-81T78?6i67^ 
-0 i?l9n?^n?3?27o ' -1'1972546577454. ' 

J;'i?5^°''30543613 , -1.70962727OfiO«ir 
-0.99022841453552 r3329787254?-?^' 

-1.5678428411484 . 0.32700437307???' 

t39l7S'f?o''|''^^^^^^^^^ 

- ^973139047623 0.66275292634S64' 

ofl3J?^^^v^^lo^^v -0.24445547163486: 

35 JJ.823533356 18973 , 0.97751 87960907 « 

2.7966260910034E-02 1 24598^780171 
0.97553628683090 . 0.74 7^5905^^^^ 
J2365598679 1.1787?07262?2V.' 

Vlt?AV.^!Sll2' 0.2«S67903l'' 
-15?02f f7??7S^ ;.^-25»^9«2231140. 
-1 493?2n^?r?? • 0-46876969933510. 

1 5572?5??olS • -0-^7981052398682 

-1.5809999704361 aieSw^^K' 
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* 1 1434341669083 , 1.5573689937592 
A 931414W81406 .0.241873 17490578. 
-1 5104018449783 . 1.2377147674561 
1 9451803353165 . 0.14569547772408. 
12988004684448 . 1.4411246776581. ^ 
w .192TO929403305 . -5.7761263102293E-02, 
-13509848117828 1.3267296552658. 
1 8577233552933 . .2.n26555278897E-02. 
-13225439786911 .-1.2418869733810. 
1:7606525421143 -9.055098 1462002E-02. 
« 1 1161555051804.-1.2882466316223. 
1 6323974132538 . -0.150467917323 11. 
0 97648859024048 .-1.2267235517502, 

0 2011785^62138 1.4773093461990. 

-0:8S89?77979660 . 1^44639 U^^^^^^^ 

1 3004480600357 . -0.24329216778278, 
^ 0 6676533^1985,-1.0373107194901. 

0 27767503261566 1.1074353456497. 
S 50618971681595 . -0.91855829954147. 
S:I§428864955902:.W^^^ 
0 79054325819016 , 0.34204551577568, 
25 -0 69724625349045 0.32767096161842, 
-0 18177013099194 , 0.65859514474869, 
0 34578496217728 0.49258133769035 

52803617715836 , 
-0 29641786217690 . 0.36104345321655, 
30 -0.40400594472885 . 0.11642142385244. 
0.11455089598894.-0.37470200657845, 
0 29129487276077,-0.24794004857540. 
57727126628160E-02 . -0.38791200518608. 
0*19418933987617 , -0.36395436525345. 
-0.40166890621185 , 0.18578046560287. 
^ -0 46209335327148 . -0.11633513122797. 
-0 41677090525627 . 0.29576727747917. 
.6!0622435063124E-02 . 0.54043716192245. 
-0 43378528952599 . 0.37475159764290. 
-0*59756731987000 , -2.909648977220 lE-02, 
-0.42078813910484 , -0.45302680134773. 
0 63260388374329 , 2.2896632552147E-02, 
0.47454643249512 . -0.43297693133354. 
0.64523047208786 . 4.222498089075 1 b-02. 
-0 41148558259010 . -0.49813264608383. 
45 8.'6346767842770E.02 . .0.63570129871368. 

0.52332174777985 . -0.35753899812698. 
-0.15361997485161 , 0.60483527183533. 
0.54942035675049 . -0.27337777614594, 
0.24155624210835 , -0.55399334430695, 
50 0.57553601264954 . -0.16218262910843. 

0.48504012823105 , 0.34690761566162. 
0.60061538219452 , -2.7974 14548695 1 £-02. 
0.46577885746956 . -0.40029266476631, 
-0.62349128723145 . -0.12451778352261. 
55 -0.30245357751846 . -0.5937591791 1530, 
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-0.64292871952057 , -0.29000273346901. 
0.72607111930847 , -0.19453555345535. 
0.65767973661423 , 0.46279S23780060. 
0.85773074626923 , -7.97747746 10996E-02. 
0.66653263568878 , 0.63700795173645. 
n <J?li^^ll!°^" • 3.8463655859232E-02, 
"^.'^J^}}^^^^^^ ^562 , -0.80670833587646. 

. -0.15678122639656. 
0.96613508462906 . -0.66218781471252. 
'^^V^^I^^P'^^^^ ' 1.1993353366852 
1.1098629236221 , -0.64718896150589. 
1.2801918983459 . 0.38049080967903. 
2-5HI!8*^3^8103 , 1.2329781055450, 
"?-«?IP^^^°3630 . 1.3380939960480. 
0.58856463432312 . 1.3312312364578. 

• -1-3701412677765. 
-1.4011698961258 , 0.54447442293167. 
i*^TlPP853851 . 0.64058619737625. 
-1.4402492046356 , 0.49066120386124. 
1.3491046428680 . 0.69823592901230, 
'2-^H^!'^31940269 .-1.4469082355499. 
-0.73903095722198 . 1.2944310903549. 
-1 .4206207990646 . 0.3558528 125286 1 
-1.2107812166214 . -0.76239341497421 

"J nJiLMI^^^^^^ ' 0.27647992968559.* 
1.0996322631836 . 0.76832592487335 
1.2723374366760 . -0.19058661162853 
' 0.75740247964859. 
n on.H^^^^^^'*^^ • -9.9505998194218E.02. 
-^'^^,^^^^^P^'2-601 , -0.730749607086 18. 
«'iiy^^'*^289825 , 4.7 180349938571 E-03 
-0.69000792503357 , 0.62838578224182 

"2-?^5?5JS^^^^O^ • -9.2191 14482402SE-02. 
0.63727533817291 .-0.43797433376312. 
'2*^?^^l°3^^'^23 . -0.18957243859768 
0.57503563165665.-0.23856720328331, 
o^El^^^^'78^506 . -0.4796301 1264801 
3.5298 123955727E-02 . 0.50607603788376 
• -0.28557923436165 • 
^ 0.16655600070953 
^^^i^^^l^.t^^'^^^^-O^ . -0.46805140376091. 
■S-??J?????7^P«^ • 0.36008602380753. 
'n??i}iiJJiZ?i?^^ • -9.266207367I818E-02. 
0.54514288902283 , -0.28926151990891 
0.26530939340591 , -0.62695705890656* 
0.71201461553574 , -0.22391571 104527 
-0.69442093372345 -0.4198 008648872 

"o'^^fllf????!^?? • -0.8580374m6660: 
-0.55144131183624 . 0.75251996517181 

0.12049946933985 0.9794833660i 257 ' 

225952148 . 0.65618216991425; 

i^-^^S}^'^^'^'^^'^^ ' -8.6988314986229E-02 
0.83770149946213 . 0.73083782196045. 
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-6.79078 17661762E-02 . -1.1372458934784, 
-0.86384069919586 . -0.77313911914825. 

1.1698449850082 , -6.432849 1687775E-02. 
-0.87872487306595 , -0.78181076049805. 
-i.1705212593079 , 7.6740458607674E-02, 
-0.75661289691925 , 0.88237941265106. 
-0.10502639412880 . -1.1395549774170, 
0.87500441074371 , 0,69834393262863, 

1.0780709981918 . -0.14845751225948. 
-0.60881644487381 , 0.85695564746857. 
0.20570844411850 , 0.98799550533295. 
-0.49079838395119 , 0.82871961593628, 
-0.87198609113693 . 0.27489036321640. 
0.34792473912239 . -0.79089194536209, 
0.35360851883888 . 0.73334228992462, 
-0.18458327651024 . 0.74415004253387. 
-0.57590019702911 . 0.43903124332428. 

5.7776970788836E.03 . -0.68923103809357. 
0.52798008918762 . 0.40390974283218. 
-0.62690776586533 . 0.18303082883358. 
0J22 19048 1424332 , -0.61702972650528, 
0.55797028541565 , -0.37609174847603. 
0.70261883735657 , 3.4565746784210E-02, 
-0J67545 17555237 , -0.48320794105530. 
-0.78116804361343, 0.15341848134995, 
-0.40339776873589 , 0.75160259008403, 
0.84919923543930 , -0.33748552203 178, 
-0.31929352879524 . 0.92272567749023. 
-0.51332873106003 . -0.90345317125320, 
0.23162299394608 . -1.0757941007614, 
-0.94100481271744 . 0.67701315879822. 
-1.2062630653381 .-0.14108607172966. 
-0.95936530828476 . 0.82507961988449. 

1.3102972507477 , 4.83580157 16076E-02. 
-0.95462840795517 . -0.95657634735107. 

1.3848868608475 , -4.5906782 150269E-02. 
-1.0633902549744 , -0.93128252029419 

1.4279329776764 , -0.14106513559818, 

1.1497660875320 , 0.88279163837433. 
0.23647691309452 . 1.4383053779602. 
-0.81110906600952 , 1.2128539085388. 
0.33149531483650 . 1.4158697128296. 
-1.2524482011795 , -0.71695524454117, 
0.42545720934868 , 1.3614822626114, 
-0.60175400972366 . 1.2690199613571, 

1.2769542932510 , -0.51767551898956, 
0.46760299801826 , -1.2636756896973, 
0.60743457078934 , 1.1649851799011, 

1.2380973100662 , 0.31721931695938, 
0.69398546218872 , 1.0290712118149, 
0.15386630594730 ,-1.1944670677185. 
-0.77655005455017 , -0.87338465452194, 
-1.8741229549050E-02. , -1.1353790760040, 
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-0.7026374936103S , 0.85432267 1S9026. 
-0.19654582440853 . -1.0637335845184, 
0.92648190259933 . 0.52192509174347. • 
0.98265725374222 , -0.37526836991310, 
-0.99220246076584 . -0.33656212687492. 
0.55053806304932 . 0.89534473419189, 
1.0506724119186 , 0.15191125869751. 
0.71802806854248 , 0.80500024557114, 
1.101 1 143922806 , -2,6790 169999003E-02. 
0.87359136343002 , 0.71470791101456, 
-0.19459256529808 , -1.1428083181381. 
-1.0133891105652 , -0.62734061479568, 
-0.34699076414108 , -1.1751155853271, 
0.54547435045242 , -1.1340159177780. 
0.48006322979927 , 1.1975069046021, 
-1 .232608 199 1196 . -0.47 1 3 1305932999. 
0.59059274196625 , 1.2095856666565, 
0.40663030743599 ,-1.3069387674332, 
0.67615872621536 , 1.2111126184464, 
1.3554904460907 , 0.35272744297981, 
0.73520946502686 , 1.2020268440247. 
-1.3775087594986 . -0.31040790677071. 
-1.1824613809586 , 0.76710045337677, 
-0.27996954321861 , 1.3730301856995, 
0.77210420370102 . 1.1527576446533. 
1.3428838253021 , 0.26121476292610, 
-0.75139236450195 , -1.1134694814682, 
1.2886705398560 , 0.25347581505775. 
-1.0653619766235 , 0.70698553323746, 
-1.2127152681351 , -0.25565749406815, 
1 .0094038248062 , -0.64 1 6779 1 604996. 
1.1179939508438 , 0.26628923416138. 
0.94675332307816 . -0.55893445014954. 
0.28359553217888 . -1.0080425739288. 
-0,46277087926865 , -0.87873333692551. 
-0.88684558868408 . -0.30556473135948 
-0.80680304765701 , 0.35761266946793, 
0.75870722532272 . 0.33003318309784, 
0.24814304709435 . 0.73252141475677. 
-0.62811541557312 . -0.35476726293564. 
0.65750586986542 . -0,13914376497269. 
-0.49959596991539 . -0.37754675745964. 
0.58338809013367 . -3.5333048552275C-02 
0.37756574153900 , 0.39624589681625. 
0.51176452636719 , 5.879380 18 143 18E-02. 
0.40890917181969 . -0.26618802547455. 
-0.44414889812469 , -0.13911458849907, 
-0.41381889581680 , 0.16923630237579 
-0.38192304968834 . -0.20201873779297 

8.9968837797642E-02 . 0.40955415368080, 
-0.24452674388885 . 0.32629245519638. 
0.39503705501556 , -3.1020458787680E-D2. 
0J26439002156258 , -0.27824369072914. 
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-0.36956810951233 , -5.6873057037592E-03, 
-0.26016449928284 . 0.23851068317890. 
-0.33284601569176 , -1.9010026007891E-02, 
-0.23125946521759 . 0.20754601061344. 
-0.28497639298439 , -8.6228372529 149E-03, 
0.17795602977276,-0.18550240993500, 
-0.22646434605122 , 2.4967700242996E-02): 
double frlt514]; . 

mainO 

15 [ 

int cn^rn:=<),crTOiim»0.nuiiibcr.mf,inii,m,id[257],s[^ 
int jrcjimj; 

double mctricp,mcuicm,miiul,mini2 jcy; 
double dcltphas, xrc, xim, cuq}lias[129]; 



20 



55 



asm(".global declared"); 
asm ("declared:"); 



nunibcr = 128; 

25 /♦♦♦«4i«#*4r4i«4i«« tijc I & Q componcnuj & find dieir ampliiudc. "frl" 

for (j=0; j <= number, j-M-) ( 

curphasIJ] = atan2(r(2*j+ll.z[2*jl); 

frl[2*j] = $qrt(z[2*j]*z[2*j]+z[2*j+l]*2[2*j+l]);} 

30 

abm(".global rcad_daia"); 
asm ("rcad.data:"); 

Esiimaic die Group delay "fr 1 " from die ampliiudc "frl " 

^ for (k=l;k<= number :k++)( 

frl[2*k-l] = (frl[2*k] - frl(2*k-2])/frlt2*k]; 
frl(2*k-2] =0.0;) i J. 

frl[2*k] = 0.0; 

^ cstimaic tlic group delay "frl" using Uic Iiilbcn iransfonn ♦••**»♦•••••/ 

m -7; 

ai-mC'.global b_ffl"); 
asm ("b_fft:"); . 
^ fft(number.ni.frl); 
asm(".global a_ffi"); 
asm ("a^fft:"): 

for(k=l; k<= numbcr/2; k++)( 

so jcy = frl(2*k-ll/numbcr, 

frU2*k-l] = 0.0-fTl(2*k]/numbcr, 
frl(2*k]=jcy;) 



for (k=numbcr/2+l; k <= number. k++) { 

1 1 
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jcy = frl[2*k-l]/numbcr; 
frl[2*k-l)= frlt2»k]/numbcn 
firl[2*k]=0.jcy;) 

asm(",globalbjm"); 
asmC'b.iffu*'); 

iffi(nunibcr»m,frl); 
asmf.global a^ifft"); 
asmrajfft:''); 

/«*****««*M«**)*4.***«i4***4i* Tlircsholding & scgnicnling ****♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦/ 

The purpose of ihis section is to deicrmine Uic segcnicnts of the received frequency frame 
whore Uic estimated group delay is larger ihan a ceriain tlireshold« Tlicse segments arc 
labelled 

id[l]=l; 
m =5 0; 
init = 0; 

for (k=l; k<»nimiben k++) { 
if(( frli2*k.ll > 0.07 frl[2*k-l] < .0.07))( 
frl[k] = frl[2*k-l]; 

if(init=0)( 

idtm] = k; 
init = I;)) 
else ( 

frl[k]=0.; 
if(init=l) ( 
m++; 
id[ml = k; 
inii = 0;)l 

if(init=l){ 
m++; 

jd[m] = number + 1 ; ) 

asmC'.global threshold**); 
asmC'ihicshold:"); 

pic purpose of tliis part of ihc software is lo dcicmiinc Uic proper sign of esiimaied group 
delay. This is done through defining two Hamming distances one if we assuem the 
cstaimied group delay was positive and llic other if it is assumed negative and choosing die 
sign that leads to the smaller disiance 

for (k=l; k<=m; k-H-){ 
mctricp = 0.; 
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mctricmssO.; 

if( ad[kl=id[k+l].l) II (id[k]=id[k+l].2) )( 
s[k)=0; 
goio(wchty;) 

for (j =id[k]+l; j<=id[k+l)-l; j-H-) ( 

dcUphas = curphas(j] - cuq>hasO-l] - &![)]; 

xrc = 1.4142l36377*cos(dcliphas); 
xim = 1.4142136377*sin(dchpha5); 

minil = 1,-xrc; 

minil =ininil*miml; 
if (mini 1>(( 1 •+xrc)*(l .+xrc))) mini 1 =( 1 ,+xrc)*(l •+xrc); 

mini2 = I. -xim; 

mini2 » mini2*mini2; 
if(mini2>((l.+xim)*(L+xtm))) mini2-(l.+xim)*(h+xim); 

mcuicp minil + miiu2; 



dcllphas = curphaslj] - curphas(j-l] + frllj]; 

xrc= 1.4142136377*cos(dcliphas); 
xim = 1.4142136377*sin(dciiphas); 

minil L-xrc; 
minil = minil * minil; 
ifCmini 1>(( 1 .-f xrc)*(l .+xrc))) mini IK 1 .+xrc)*(l .+xrc); 
mini2 = L-xim; 
mini2 = mini2 * mini2; 

if(mini2>((l.+xim)*(l.+xim))) mini2=(l .+xim)*(L+xim); 

mctricm+= minil +mini2; 

} 

5[k]=0.1; 

if(mctricp < mclricm) 5fk]=l ; 
/♦ prinifrs[%d] = %dVn",k,s[k]);V 

twenty: k++;) 

for (k=l ; k<=m; k+H-)( 

for 0 =id[k]; j<=id[k+lM; j++) frlfjj s[k]; 
/*k++;*/ ) 

asm('\gIobal a_f;>"); 

asm("a,fs:"); 

Decode the Signal 
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15 



20 



25 



30 



1 = 1; 

while 0 < numbcr+l){ 

dcllphas = curphas[I) • cut)1us[H] - frl[l]; 

xrc =s cos(dcltphas); 
xim — sin(dcUphas); 

jrc = xrc/fabs(xn:); 
jim = xini/fab5(xim); 
if (jrc = jim) ( jrc = 0 -jrc; 
jim = 0 -jim;) 

if(I>l)( 
ifOrc != Jj2[2*l]) cm>rrc++: 

ifQim != jj2[2*l+l]) cnorim++; 

) 

printf('*%d %d %dSii",i,crrorrc.crrorim); 

asmCglobal end"); 
asmC'end:"); 

) /*cndof main*/ 



35 



Claims 



1- A transceiver including a transmitter for transmitting electromagnetic signals and a receiver for receiving 
electromagnetic signals, the transmitter comprising: 

^ an encoder for encoding Information; 

a wideband frequency division multiplexer for multiplexing the information onto wideband frequency 
channels; 

a low pass filter; 

a local oscillator for upconverting the multiplexed Information for transmission; and 
45 a processor for applying a Fourier transform to the multiplexed Information to bring the information 

into the time domain for transmission. 

2- The transceiver of claim 1 in which the received signals have amplitude and phase characteristics and 
further including, in the receiver of the transceiver 

50 a channel estimator for estimating one or both of the amplitude and the phase differential of the 

received signals to produce as output one or both of an estimated amplitude and an estimated phase dif- 
ferential respectively. 

3. The transceiver of claim 1 or 2 further including, in the receiver of the transceiver 
55 a bandpass filter for filtering the received electromagnetic signals; 

a local oscillator for downconverting the received electromagnetic signals to produce output; 
a sampler for sampling the output of the local oscillator to produce sampled signals to the channel 
estimator; and 
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a decoder for producing signals from the sampled signals and the output from the channel estima- 
tor. 

The transceiver of daim 3 further including, in the receiver of the transceiver: 

a deprocessor for applying an inverse Fourier transform to the samples output from the sampler. 

The transceiver of claims 3 or 4, further including, in the receiver of the transceiver. 

a power controller before the sampler for monitoring the power of the received signal and for con- 
trolling the power of the signal. 

The transceiver of any one of claims 2-5, further including, in the transmitter of the transceiver. 

a pre-distorter before the processor, the pre-distorter being connected to the channel estimator, 
for predistorting a signal to be transmitted with one or both of the estimated amplitude or the estimated 
phase differential. 

J. The transceiver of daim 6 in which the power controller is also connected to the pre-distorter for controlling 
the power of the signal to be transmitted. 

8. The transceiver of any one of daims 2-7 further induding: 

means to modify the received signal with one or both of the estimated amplitude and phase dif- 
20 ferential respectively. 

.9. A method for allowing a number of wireless transceivers to exchange frames of information, the method 
comprising the steps of: 

multiplexing a first frame of information over a number of frequencies within a frequency band at 
25 a first transceiver to produce multiplexed information; 

processing the multiplexed information at the first transceiver to produce processed information; 
transmitting the processed information to a second transceiver using a carrier frequency fc; 
receiving the processed information at the second transceiver; and 

processing the processed information at the second transceiver during a first time interval; 
30 in which the width of the frequency band is chosen so that neither carrier nor clock recovery is 

required at the second transceiver. 

10. The method of daim 9 in which the frequency band is formed from a first set of Ki points and a pair of 
tail slots each having K2 points, each of the points being separated by a frequency range Af, the second 
35 transceiver has a maximum expected clock error xT, vyhere T is the duration of one time domain sample, 

the information is multiplexed over a number M of levels, and is selected such that 29cx/Ki ^ n/M, where- 
by the width of the frequency band is chosen so that neither carrier nor clock recovery is required at the 
second transceiver. 

^ 11. The method of daims 9 or 10 further including transmitting a second frame of information from the second 
fransceiver to the first transceiver within the same frequency band. 

12. The method of claim 11 in which K2 is selected so that the out of band signal is less than a given level. 

13. The method of claims 11 or 12 in which the first and second transceivers have an expected maximum 
^ relative velocity V, the first and second transceivers have carrier frequencies with a frequency offset froni 

each other of fo, the carrier frequency has a corresponding travelling wavelength X and Af is selected so 
that [V/(XAf) + fo/Af ] is less than or equal to a preselected sampling error. 

14. The method of any one of claims 9 - 13 in which processing the rnultiplexed information at the second 
^0 transceiver further includes calculating the mean of the phase shift due to sampling error by summing 

the estimated phase differential of the received signal. 

15. The method of claim 14 in which the mean of the phase shift due to sampling error is divided by Ki and 
the result removed from the phase differential of the received signal. 

55 . ' 

16. The method of any of claims 11-15 further comprising: 

transmitting the second frame of information from the first transceiver to the second transceiver 
after a preselected time interval after transmitting the first frame of information from the first transceiver 
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to the siecond transceiver; and 

during the preselected tinne interval, exchanging information between the second transceiver and 
at least one other transceiver. 
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